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SUMMARY 


A two-dimensional finite-difference numerical technique 
is presented to determine the temperature distribution in a 
solid blade of a radial turbine guide vane. 

A computer program is written in Fortran IV for IBM 
370/165 computer. The computer results obtained from these 
programs have a similar behavior and trend as' these obtained 
by experimental results. 
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INTRODUCTION 


The power output of a radial gas turbine is directly pro- 
portional to the inlet temperature of the combustion gases. 
Generally, higher inlet temperatures are desired in order to 
achieve higher power outputs . Due to the limitations on the 
material strength, the operating temperature is limited for 
the safe operation of the blading without failure r One way 
of ensuring the safe operation of gas turbine blades is to 
provide adequate cooling for the blades. The cooling of the 
blades helps the safe operation of an engine by reducing the 
temperature level in the blade material and by equalizing the 
temperature differences throughout the blade sections. By 
determining the temperature distribution in the blades of a 
gas turbine, the critical stress areas could be located and 
hence could be provided with sufficient cooling mechanisms. 
Kuhl^^^ reported the temperature measurements that were taken 
in the rotor blades of a gas turbine. He used the analogy 
between the heat flux and the electric current within a 
three-dimensional model to determine the heat flow in a 
complete blade. 

Calculating the temperature distribution at a cross section 
of a gas turbine blading using analytical methods poses some 
problems due to the irregular shape of the boundary of the 
blade and the fact that the temperature and the velocity of 
the gas around the blade surface are variable from one point 
to another. The variable gas velocity around the blade 
results in a variable heat transfer coefficient. Due to 
these circumstances, it is necessary to resort to numerical 
techniques using a finite-difference or a finite element 
method of solution. In the present study, a computer program 
is Written using a two-dimensional finite-difference numerical 
technique to determine the temperature distribution at a 
cross section of a solid blade. Data concerning the geometry 
of the blade surface is given as an input in the form of the 
coordinates of the points intersection of the mesh lines of 
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Figure 1 with the boundary of the blade. The gas temperatures 
and the corresponding heat transfer coefficients at these 
blade boundary points are also given in the input data. 

This computer program determines the two-dimensional 
temperature distribution in the cross section of a blade without 
cavities. It can be modified however to handle the case of 
a hollow blade utilizing internal cooling as in Reference [2] . 

It is assumed that the thermal conductivity of the blade 
material is the same everywhere since the variations in the 
conductivity will not be significant except for very large 
temperature differences within the blade. A subroutine can 
be added to the program however if it is desired to take 
into account the local variation of the thermal conductivity 
with temperature. The second order derivatives in the partial 
differential equation governing the temperature distribution 
in the blade are expressed in terms of second order accurate 
finite-differences and the resulting algebraic equation in the 
mesh point temperatures are solved using the Gauss-Seidel 
iteration technique. Over-relaxation factors are used to 
accelerate the convergence process. 


GOVERNING EQUATIONS 


In plane rectangular coordinate system, the differential 
equation governing the temperature distribution in the turbine 
blade is given by: 



( 1 ) 


In the above equation it is assumed that there is no heat 
generation in the blade body and that the thermal conductivity 
of the material is constant. 

The boundary condition at any point on the blade surface 
can be expressed as; 


K = h(T - T„) 
0 IX ^ S 


(2) 
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Where 


K 


is the thermal conductivity of the blade material 

n is the normal to the blade surface at a given point 

Tg is the gas temperature at the surface of the blade 

and Tg is the blade surface temperature at the particular 
point in question 

The turbine blade and the grid network used in the present 
case is shown in Figure 1. The particulars of the blade are 
given in Table 1. The grid spacings in the x-direction, DX, 
are uniform as well as the spacing DY in the y-direction. 

Referring to Figure 2, the partial differential Equation (1) 
can be written in a finite-difference form at any mesh point 
P(i,j) inside the blade surface in terms of the temperatures 
at the neighboring points A, B, C and D as follows: 


T. T - 2T. . + T.^t . T. . T - 2T. . T. 

1-1,3 J-O i+lQ + ^>3 + ^ = 0 

2 2 
DX^ DY 


or 


T. 


2(DX^ + DY"^) 


(3) 


Equation (3) is second order accurate, and is valid at all 
interior points whose four closest neighboring points happen 
to be at regular mesh locations. For points near the boundary 
of the blade, in Figure 3, the neighboring points A, B, C and D 
do not occupy normal mesh locations and hence the expression 
given by Equation (3) has to be modified to take into account 
the nonsymmetrical locations of the neighboring mesh points. 

In a general case such as in Figure 4, where points A and C 
are located at and r;^ fractions of DX and points B and D are 
at 6 ^ and 62 fractions of DY relative to P, the difference 
Equation (4) will be used. This equation is arrived at by 
expressing the temperatures at A, B, C and D in terms of Taylor 
series expansions at P, and then eliminating the first order 
derivatives between these expressions to obtain the following 
equation: 


A 



DX^(^^+?2) ^ 


- T, ^ (~ + ^) + — 


ifD 


D ‘^( 6 j^+ 62 ) 


[!i|l_i - T. . (^ + ^) + = 0 (4) 

Ot ^ f 3 °o °o 


Equation (4) can be written in a more convenient form as 
follows: 


3-fD 


1 

E 


T, 


T. 


T. 


T. 


r i-lf j . "if j+i ^ /DX^ ^ i/j-1 ^ _ i fj+1 T 1 


Where 


(5) 


E = [ 


4^2 


+ (^) 


'l^2- 


( 6 ) 


The above equations are second order accurate when “ '^2 

and ~ ^2 * equally spaced grid lines, the factors 

Cl, C 2 f and 62 >^ill all be equal to 1.0 and Equations 
(5) and ( 6 ) will ultimately reduce to Equation (3) . Equation 
(5) is used instead of Equation (3) at grid locations next to 
the boundary of the blade such as at the point P in Figure 3. 
At the blade boundary points such as C or D in Figure 3, the 
convective boundary condition given by Equation (2) has to be 
used to evaluate the surface temperature of the blade. In 
order to apply Equation (2) , the normal gradient of the 
temperature normal to the boundary point has to be expressed 
in a finite difference form. Referring to Figure 3, if NON' 
is the normal to the blade surface at the boundary point D, 
then the temperature gradient D is approximated by: 


T - T h 

9T ^ ^ ^ - 3 fT - T ^ 
9n DN " K ^ g 
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Rearranging Equation (7) , the temperature at the boundary point 
D is given by: 


h DN 

[T + — T ] 

D DN 


[1 + 


‘D 


K 


-1 


( 8 ) 


The solution of Equation (3) at the interior points, (5) 
at the points next to the boundary, along with Equation (8) 
at the boundary points will give the temperature distribution 
in the blade. 

COMPUTER PROGRAM 

The details of the input to the computer program are given 
in Appendix A. Prom the blade boundary points, the program 
calculates data concerning the 

mesh structure and the boundary point coordinates. As an 
initial guess the temperature matrix T is set to any convenient 
value. The Gauss-Seidel method of iteration then starts with 
the interior points of the blade. The computations giving new 
temperatures of mesh points on each I line proceed in the 
X direction. Equation (3) is used for all points which do not 
lie next to a boundary and Equation (5) is used for all the 
points which lie next to a boundary. Once this is completed, 
the temperatures on the blade surface are evaluated using the 
convective boundary condition expressed by the finite-difference 
Equation (8) . The normal to the boundary at any point is 
determined by fitting a least squares parabolic curve through 
the three boundary points consisting of the one under considera- 
tion and the two adjacent points on both its sides. Prom the 
parabolic equation of the curve, the slope of the tangent to 
the surface at the point can be found and hence the slope of 
the normal to the curve at the required. The coordinates 
of the intersection of the normal with the closest mesh line 
inside the blade, point N in Figure 3, can thus be found. The 
temperature at the intersection point, N, is determined by 
linear interpolation from the temperatures of the neighboring 
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points A and P of Figure 3. The length DN is determined from 
the coordinates of N and D which are known from previous 
steps. In each iteration, the boundary temperatures are first 
evaluated at points on the I lines and then the temperatures 
at the boundary points on the J lines are determined in a 
similar fashion. This iteration process is repeated until 
the sum of the squares of the differences between two successive 
iterations is less than a prescribed small quantity. 


i.e. I (T - T <_ e C9) 

all points 

th 

Where T^ is the temperature after n^ iteration 

T^_^ is the temperature after (n-l) th iteration 

and e is the prescribed error limit. 


The process of convergence could be accelerated by using an 
over-relaxation factor, oj, thus: 

T = 0) t' + (1 - w)T„ , (10) 

n n n— i 


Where is the temperature after (n-l) iterations 

I 


T^ is the temperature computed after n iterations 
(ii is the over-relaxation factor 

5 1 

T^ is the temperature used for the (n+1) iteration. 

After the final iteration, all the grid temperatures are printed 
and the locations for isothermals are found by linear interpolation. 


COMPUTATION PROCEDURE AND DISCUSSION 

To use the program for determining the temperature dis- 
tribution, the input data has to be prepared. The details of 
input parameters required are given in Appendix A. Defining 
the blade boundaries can be most easily done graphically as 
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shown in Figure 1. To find the gas temperatures and heat 
transfer coefficients around the blade, the flow behavior 
around the blade should be obtained first. 

The velocities on a radial gas turbine blade are deter- 

(3) 

mined from the computer program given by Dastanis . This 
can be used to determine the flow properties on a blade-to- 
blade surface in a turbomachine. The solution is obtained using 
a velocity gradient method; which uses information obtained 
from a finite-difference stream function solution at a reduced 
weight flow. From the velocity distribution, the temperature 
of gas around the blade can be determined. A typical velocity 
distribution for the case used in this investigation is shown 
in Figure 5. For convenience, the blade surface flow 
velocities are plotted as a function of surface distance; which 
are measured from the point, in Figure 1. The temperature 
distribution is computed using the energy equation and is 
shown in Figure 6 as a function of the blade surface distance 
from the leading edge. The heat transfer coefficient is 
determined from the Reynold,' s Analogy and the boundary layer 
characteristics . 

The boundary layer solution is determined using the computer 

(4) 

program written by Herring and Mellor' ' . The velocity 
distribution of Figure 5 is used as input for the boundary 
layer computations. The details of the program input and 
method of solution are given in Reference (4) . This program 
can be used for both laminar and turbulent boundary layers 
over two-dimensional or axi-symmetric bodies. The program 
also takes into consideration the effects of pressure gradients, 
Reynolds Number, wall transpiration and surface roughness. 

This program solves the boundary layer equations numerically 
by using Crank -Nicolson method in combination with a fourth 
order Runge-Kutta solution scheme. In its output, the program 
gives the boundary layer parameters on the blade surface, such 
as the displacement thickness, the momentum thickness, the 
friction coefficient, and the energy thickness. Knowing 
the friction coefficient along the surface from the boundary 
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layer computations, the heat transfer coefficient is determined 
using Reynold's Analogy. The convection heat transfer 
coefficient shown in Figure 7 is calculated using the gas 
velocity and temperature distributions of Figures 5 and 6. The 
temperature and the heat transfer coefficients around the 
stagnation point are not expected to be accurate, however, 
since the velocity distribution calculated using the method 
of Reference (3) is not accurate in this region. 

The present conduction heat transfer program uses a 
40 X 12 mesh point with DX = 1.27 mm and DY = 0.635 rrmi. 

The over-relaxation factor for faster convergence 
is determined by trial and error. Figure 8 gives the 
variation of the number of iterations needed for a certain 
convergence with the relaxation factor, w. Three different 
error limits, e = 0.1, 0.05, and 0.02 were considered for 
convergence. With an error of e £ 1, the minimum number of 
iterations occurs at a> = 1.75, Similar behavior was 
observed with the other error limits, with the minimum number 
of iterations being obtained with to = 1.7 and 1.6 for 
e ^ 0.05 and 0.02 respectively. The solution was found to 
oscillate at higher values of the over relaxation factor to. 

The program written in FORTRAN IV for an IBM 370/165 
computer, is given in Appendix B. The computation time was 
14 cpu seconds for 280 iterations. An explanation of the 
program input and its preparation is given in Appendix A 
while in Appendix C, a sample output is given. Figure 9 
shows some of isothermals for the blade profile considered. 

From Figures 6 and 9 it can be observed that the upper surface 
of the blade is generally at a higher temperature than the 
lower surface. The isothermal lines, which are obtained 
using the present numerical methods are found to have similar 
behavior and trends as those found experimentally by- Kuhl 

CONCLUSIONS 

The computer program presented here gives a two-dimensional 
theoretical temperature distribution for a radial turbine guide 
vane blade. The temperature distribution agrees with the 
experimental results obtained by other investigators. 
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LIST OF SYMBOLS 


Points indicated in Figure 2 

Distance along normal to surface of blade (Figure 3) 


DX 

DY 

E 

h 

I 

J 


K 


Grid spacing in x direction (m) 

Grid spacing in y direction (m) 

Constant defined by Equation ( 6 ) 

Heat transfer coefficient C 5 ) 

hr 

Grid line numbers in x direction 
Grid line numbers in y direction 

Thermal conductivity of material of the blade ( ) 

hr m ®K 


n 


P 

T 


X 

y 



^2 

e 

Subscripts ; 


Normal direction to the blade surface 

A point of grid structure as in Figure 2 

Temperature 

X-coordinate 

y-coordinate 

Fraction of grid spacing for point B in Figure 4 

Fraction of grid spacing for point D in Figure 4 

Fraction of grid spacing for point A in Figure 4 

Fraction of grid spacing for point C in Figure 4 

Ejrror term as defined in Equation C9) 


D Corresponds to boundary point in Figure 3 
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Corresponds to gas 

Mesh line numbers in x and y directions 

Corresponds to the point where normal NN^ to the 
boundary at D cuts the mesh line AP (Figure 3) 


Surface 



TABLE I 


PARTICULARS OF RADIAL TURBINE NOZZLE BLADE: 


Leading edge Radius = 2»07 rnm 

Trailing edge Radius = 0.40 mm 


c = 49 . 3 mm 


5 

Thermal conductivity = 0.7476 x 10 ^ 

Boundary points : 


X mm 

y (Lower) mm 

y (Upper 

0.0 

2.54 

2.54 

2.54 

0.71 

4.57 

5.08 

0.63 

5.16 

10 . 16 

1.73 

5.92 

15.24 

2.67 

6.32 

20.32 

3.38 

6.48 

25. 40 

3.76 

6.35 

30.48 

3.94 

6.04 

35.56 

3.73 

5.46 

40.64 

3.22 

4.65 

45.72 

2.54 

3.63 

49.27 

2.54 

2.54 


DX = 1.27 mm ; DY = 0.635 mm. 


J 

m °K 
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1 ^ 


H 

ee. =X3(6) SURFACE DIST. OF So FROM NN =1.016 mm 

ef. =X4(6) 

cf =X5(6) DX = 127 mm 

eg =X6(6) DY =0.635 mm 


FIGURE 1. TURBINE BLADE WITH MESH LINES 



FIGURE 2 GRID STRUCTURE FOR 
INTERIOR POINTS 
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FIGURE 3. POINTS NEAR A BOUNDARY 


16 




FIGURE 4 GENERAL GRID POINT WITH 

UNEQUAL SPACINGS 
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FIGURE 5. VELOCITY DIST. AROUND THE BLADE 




FIGURE 6 TEMPERATURE DIST. AROUND THE BLADE 




FIGURE 7 HEAT TRANSFER COEFF. AROUND THE BLADE 



OF ITERATIONS 




22 


FIGURE 9 ISOTHERMALS IN THE BLADE 






APPENDIX A 


Input for the computer Program: 

The input to the computer program consists of the following 
parts : 

1. Physical parameters for the grid system/ blade 
material properties, allowable error, and number 
of iterations. 

2. Blade geometry. 

3. Convection heat transfer coefficients around the 
blade at the mesh line intersections. 

4. The gas temperatures around the blade at the 
mesh line intersections. 

The first part consists of the following items: DX, DY, XK, 

0ME, NX, NY, SUMM, NITE, NTE, N4 , N5 , NN, Nl , N2, N6 . These items 
are specified according to the Format 4F5.3, 212, F6.3, 213, 612. 

DX: The grid spacing in the x-direction (m) . 

DY: The grid spacing in the y-direction (m) . 

XK; Thermal conductivity of the material of the blade ) » 

0ME: Relaxation Factor (If not known give a value of 0.0) , 

and the program would assume a value of 1.0) . 

NX: The number of grid lines in the x-direction 40). 

NY; The number of grid lines in y-direction 20) . 

SUMM; Maximum allowable error, e, which is the sum of the squares 

of temperature differences between two successive iterations. 

NITE: Maximum number of iterations. This would terminate the 

computations if the required convergence is not achieved 
at or below the specified number of iterations. 

NTE: Print out of temperature required at every NTE iteration. 

N4: The first J line that intersects the blade in more than 

two points. Specify N4 = 99 if the above phenomenon 
does not occur. (In the example of Figure 1, N4 = 5) . 


23 


N5: Final J line that intersects the blade in four points. 

When N4 = 99, N5 represents the J line which differentiates 
the top and bottom portions of the blade near the trailing 
edge of the blade. In the example of Figure 1, N5 = 7. 

NN; This corresponds to the J mesh line which differentiates 
the top and bottom portions of the blade near the front 
end of the blade. In Figure 1, NN = 5. 

Nl: This corresponds to the I mesh line, which separates the 

two points 'a' and 'b' on the bottom blade surface inter- 
sected by the J line: N5. (If N4 = 99, take Nl = NX) . 

In Figure 1, Nl can be any number between 22 and 27. 

N2: This corresponds to the I mesh line upstream of the point 

C in Figure 1, which is the last point on the top surface 
intersected by the N5 mesh line. If N4 = 99, N2 will be 
equal to NX. 

N6: Corresponds to the J value at the point d of Figure 1, 

which is the center of the trailing edge radius . 

The second part consists of three items. The first set of 
cards contain the y coordinates where each I mesh line intersects 
the blade. The first coordinate, Yl(l), refers to the lower 
surface of the blade and Y2(I) refers to the upper surface of 
the blade, both are given in inches according to the Format 
12F6.3. These values are given in pairs of coordinates corre- 
sponding to the bottom and top boundaries at each I mesh line. 

If the I line does not intersect the blade, give the two 
coordinates as (NY-l)DY. The program reads in NX pairs of 
nxambers with 6 pairs to a card, with a total cards read = NX/6 
or the next whole number if NX is not of multiple of six. 

The next set of cards contain the x-coordinates where each 
J grid line intersects the blade surface where X3(.J) refers 
to the point closer to the Y axis and X4(J) refers to the 
next intersection further from the y axis. Both are given in 
inches according to the same Format, the 12F6.3 as before.. 

The program reads in NY pairs of numbers with 6 pairs to a card 
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with a total of NY/6 cards read. If N4 is not equal to 99, 
the third set of cards are given, which specify the other pair 
of x-coordinates for the mesh lines J = N4 to N5, These x- 
coordinates are given as X5(J) and X6(J) in inches and according 
to the same Format 12F6.3. If N4 = 99, the third set of cards 
does not exist. The various parameters are shown in Figure 1. 

The third part of the input is the convection heat transfer 

2 

coefficients in J/hr m ®K . These values are given in three 
sets of cards. The Format used here is 12F6.3. The first set 
of cards contain all H(I) values (for 1=1, NX) corresponding 
to the top surface. The computer program stores the bottom 
surface coefficients in locations H(l,l) and top surface 
coefficients in the locations H(I,2). The second set of cards 
corresponds to all HX(J) values (for J=l, NY) corresponding to 
to the nearer surface and then all the HX(J) values corresponding 
to X4(J) points. These values are given in the same units and 
according to the same Format as before. The computer program 
stores these coefficients in the locations HX(1,J) and HX(2,J) . 

If N4 is 99, the third set of cards are omitted. Otherwise 
the third set of cards are given corresponding to HXl(J) values 
and HXl ( J) values for J = N4 to N5. The HXL(J) corresponds 
to X5(J) point and HXl(J) to X6 ( J) . 

The fourth part of the program input gives the gas tempera- 
ture values in ®K. These values are given in three sets of 
cards similar to part three. These are given according to the 
Format 10F8.2 with 10 values to a card. The first set of cards 
contain all TG(I) values corresponding to bottom surface and 
then all TG(I) values. These values are stored in the computer 
in locations TG(I,1) and TG(I,2). These two correspond to 
H(I,1) and H(I,2) respectively. The second set contain all 
TGX(1,J) values and then all TGX(2,J) values corresponding to 
X3(J) and X4(J) points respectively. Finally, if N4 is not 
equal to 99, the third set of cards are given containing TGX1(.1,J) 
values for J = N4 to N5. This set is omitted if N4 = 99. 
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APPENDIX B 


PROGRAM LISTING 

In the following pages 27 to 39 is a listing of the 
main program and the subroutines SLO and BUN. 
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IV G LEVEL 21 


MA [N 


DATE = 7A273 


15/22/49 



01 MENS ION T {'+0,20) , TX( 2, ?0) , TXU 2,20) ,H( 40,2) ,HX{ 2 , 20) ,HXl (2 ,2 0) ,T 
1G{ 40,2) tlGXI 2,20) ,TCX1( 2,20) , SI (40f 2 )»S2 (40,2), Dl(40»2),D2( 40,2) ,Y 
21(40),Y2(40),X3(20) ,X4(20) ,X5(20) ,X6(20) ,IN(2,40), IM(2,20) ,IMlt2,2 
30) 

DIMENSION SK20) 

CALL UNOFLw 
WRITE(6,45l) 

451 FORMAT!//, 40X, 'TEMPERATURE DISTRIBUTION IN A TURBINE BLADE',///) 

RE AO (5, 1) C'X,DY, XK,OME ,NX,NY,SUMM,NITE,NTE,N4,N-5,NN,N1,N2,N6 
1 FORMAT (4 F5. 3, 212, F6. 3, 21 3,512) 

C READS IN INPUT DATA BOUD. PTS., ETC. 

IFIOME.EQ.O. ) QME=1.0 

IF(NX. GT.40.gr, NY. GT. 20) GO TO 505 

VlRITE(6,3) DX,DY,OME,XK,NX,NY,SUMM 

3 FORMAT! // ,5X,' OX = • , F8.4, 5X, 'DY = • , F 8. 4, 4X , * OMEGA =• ,F6.3, 4X, * THER 

IM. COND. =*,F8.4/,5X,*NX = • , I 8 , 5X , • N Y = ' , I 8 , 4X , * MA X . 6RR,= • , F8 .4 / ) 
NXl=NX-l 
NYl=NY-l 

READ(5,5) lYK n ,Y2(I) ,I=1,NX) 

5 FORMAT! (12F-6.3) ) 

READ{5,5) (X3IJ) ,X4( J) , J=1 ,NY) 

IF(N4.EQ.09) go TO 26 
READ(5,5) ( X5{ J) , X6{ J) , J=N4,N5) 

26 CONTINUE 

DO 6 1=1, NXl 
XM=(Y1( I) +0.0001 )/0Y 
LM=XM 
S 1 < I , 1 ) = 1 . 

S2( I,l)=l. 

S1U,2)=1. 

S2(T,2)=1. 

D2(I,1)=1. 

D1 (1 ,1 )=FICAT(LM+1 )-Yl ( I )/DY 
IN(1,I )=LM+2 
XM={Y2(I)+O.OOOl)/DY 
LM = XM 

OK 1,2 )=1. 

IN(2,1 )=LM+1 

IFCAB5(Y2( I)-FLQAT(LM)«DY).LT.0.0001) IN(2, I) = LM 
D2( I,2)=Y2( I )/OY-FLCAT( INI 2, I)-l) 

6 CONTINUE 

30 FORMAT(5X,I3,2( I5,F10.A) ) 

DO 7 J=2,NYI 
XM=(X3I J)+0.0001)/DX 


LM = XM 

IM(1,J)=LM+2 
JJ=IM( 1 , J) 

si( JJ,i )=float(lm+i)-x3( j)/nx 
IF< J.GT.NN) Sl( JJ,2)=S1( JJ,1) 
TF(J.GT.NN) SKJJ,l)=l. 

XM = (X4( J) +0.0001) /DX 
. LM=XM 

IM(2,J)=LM+1 
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IF {ABS(X4( j) -FLOAT (LM)»!'DX)*LT.0.0001) IMI2f J)=LM 
JJJ=IM(2fJ) 

S2 { JJ J , I ) =X4 (J ) /DX-FLO AT ( J J J- 1 ) 

IF(J,GT.N5) S2( JJJ »2)=S2( JJJ, 1) 

IF(J.GT.N5) S2(JJJ,1)=1. 

7 CONTINUE 

35 F0RMATI5X,I5,2( I5,2F10.4n 
Sl(2fl)=l.-X3(IN(li2))/OX 
IA=IN(l,2)+l 
IB=IN(2,2 )-l 
DO 9 I = IA»Iti 

9 S I( I-IA + 1 )=1 .-X3n A)/DX 
1F(N4*E0.99) GO TO 12 
DO 8 J=N4fN5 

XM= I X5(J) +0.00001) /OX 
LM=XM 

IMKl, J) = l.M+2 
JJ=IMU l» J) 

Sl( JJtl ) = F'LQAT{LM+ I )-X5( J)/OX 

XM=(X6(J) +0.00001 }/0X 

LM=XM 

IMl(2f J ) = LM+1 

IF(X6< J) . EO.FLOATI LM)«DX ) IM1I2, J)=LM 
JJJ=IM1(2, J) 

S2UJJt2) = X6( J )/DX-F|.nAT ( JJJ-1) 

8 CONTINUE- 
12 CONTINUE 

10 FaRMATl40l2) 

15 FORMAT! 16F5. 3) 

RE AO! 5, 20) ! !H! I,J),l=l,NX), J=l,2) 

20 FORMAT!12F6.3) 

READ! 5, 20) ! !HX! 1 1 J) fJ = l ,NY) ,I = 1»2) 

IFIN4.E0.99) GO TO 33 

READ!5f20) ( !hX 1! I,J)» J=N4tN5),I=lf2) 

33 CONTINUE 

READ! 5, 25 ) {!TG!IfJ)tI=ltNX),J=l»2) 

25 FORMAT! rOF8. 2) 

READ! 5,25) ( ! TGX !I , J ) , J= 1 ,NY ) tl = l,2) 

IF!N4,EG.99) GO TO 34 

RE AD! 5,25) ( !TGX1! ItJ) , J = N4t N5) t I = lt 2) 

34 CONTINUE 

C SETS INITIAL VALUES AND 30UNDARY POINTS 
OXY=DX*DX/Oy/DY 
00 100 I=2,MX1 
IA=IN!1,I)-1 
T3=IN!2»I)+1 
DO 100 II=lAfI0 
100 T! I f I I ) = 1900. 

DO 110 J=2,NY1 
TX! ltJ)=1900. 
no TX! 2,J)=1900. 

IF!N4.EQ.99) go TO 36 
DO 115 J=N4,N5 
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TXU It J) = 1900. 

115 TXK2, J)=1900, 

36 CONTINUE 

DO 120 1=1, NX 
IA=IN( Itl )-2 
IB=IN(2,I)+2 
DO 120 J=1,MY 

IF( J.GT.IA.AND.J^LT.IB) GO TO 120 
IF(J.LE.IA) ni,J) = TG(I,l) 

IF(J.GE.IB) T(l , J )=TG( I .2) 

120 CONTINUE 

WR1TE(6,456) 

458 FORMAT(lHl) 

WRITE(6,455) 

455 FOrMAT(/,5X, 'TEMPERATURE AT GRID POI NT S • , // , AX , * J= * , 5X, * I ' , 9X, * 2 • , 

l9X,»3»t9X,*4't9Xt*5',9X,*6',9X,'7*,9X,*8*,8X,*9»,8X,'10*t8X,*ll*,8 
2X, » l2'/t4Xt* I*/) 

DO 124 1=1, NX 

124 WRITE(6t440l I , ( T ( I , J) , J = 1 ,NY ) 

c iterations... lower, middle and top points 

WRITE(6,456) 

456 FORMaT(//, 5X,MB2,C2) IS THE POINT ON THE BCUNDARY AND (Bl,Ci),(B3 
i,C3) ARE the Surrounding points on BLA0E*t/,5x,» a,b,c are the coe 
2FFICIENTS OF THE PARABOLIC CURVE Y = A*B>!‘X+C*X’!‘X * / ) 

SUM=0. 

iter=i 

125 CONTINUE 
AL=H(l,l)/XK 

T( l,5) = (DX’!‘AL*TG( 1,1) + T(2,5) )/(l.+AL*DX) 

DO 150 I=2,NX1 
IA=IN(1,I)+1 
IB=IN(2,I }-l 

A=l./(Sl{ T,1)«S2 Ut 1) )+4. / (OKI, 1)«02( I,in 
Al=T{I-l, 1A-1)/S1( 1,1) /( 51(1,1) +S2( 1,1)) 

A2=T( I + l, IA-l)/S2( 1 ,1) /{ SI ( r , 1)+S2( I , 1) ) 

IFII.GT.Nl) GO TO 130 

IF(Sl(I ,1 ).M6.l.) A1=TX(1, IA-1)/S1(I ,1)/(S1( 1 , 1 ) +S2 ( I , 1 ) ) 

IFtS2( 1 , 1 ) .NE.l. ) A2=TX(2, I A- 1 ) / S 2 ( I , 1 ) / ( SI ( 1 , 1 ) +S2 ( 1 , 1 ) ) 

GO TO 135 
130 CONTINUE 

IF(S1( I f 1 ) .NE.l. ) A1 = TX1(1,TA-1)/S1( I,U/(SK I,l)>S2(i,l)) 

IF(S2( 1,1 ). NE.l.) A2=TXl (2, lA-l )/S2( I, i)/(Sl( I, 1) + S2( I, 1) ) 

135 A3=T( I , IA-2)/ni( I , 1 )/<0i( I, I )+D2( I, U ) 

A4=T ( I , lA) /D2( I,l )/(ni( I , I ) + D2( I ,1 ) ) 

TNEW=( A1+A2+0XY*( A3+A4) ) /A 

TNEW=( l.-OME)«T( I , IA-1 ) + ave.-|.TNEW 

5UH = SUM+(T( I ,I A-1 )-TNEW)»MT( I, IA-1)-TNEW) 

T( I,IA-1 j=TNEW 
IF(IA.GT.IB) GO TO 142 
DO 140 J=IA,IB 
IF(I.EQ.2) on TO 138 

TN6W=0.5<=(T( I-1,JM-T( 1+1,J) + dXY>MT(I,J- 1)+T{ ) )/(l.+DXY) 

TNEW=( l.-UHE)^T( I, J) + OME'!'TNEW 
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SUM=SUM+m I , J)-TN6W)* tn I,J)-TNEWl 
Tl I,J)=TNEW 
GO TO 140 
138 CONTINUE 

S=SU J- lAf 1) 

TNEW=( l.-OME)*T( I » J )+OME^TNEW 

TNEW=S*(TXU ,J)/( S*( l.+Sn+T ( I ♦! , J ) / ( 1 . + S I +0 .S^^-OXY# ( T ( I , J- 1 ) + T II t J 
1+1) n/( I. + S I'OXY) 

SUM=SUM+(T(I , J)^TMEW)*IT ( I ,J)-TNEW) 

T(I,J)=TNEW 
140 continue 

142 A=l./< SU I, 2)*S2( I,2))+DXY/(Din , 2)^‘02C 1,2) ) 

A1 = T( I-l, IB + lJ/SlU ,2) / (SU 1 ,2)+S2(I,2 3 ) 

A2 = T( U1,1B+U/S2(I,2)/CS1U,2)+S2(I,2)) 

IF( I,GE.N2) GO TO 145 

IF( Sl( I ,2>.NE.l. ) A1=TX( 1, IB+1)/S1(I ,2)/ (SI ( I,2)+S2 ( 1,2) ) 

IFIS2( I ,2 ) ,NE.l. ) A2=TX(2, IB + 1)/S2( I ,2)/(Sl( I,2)+S2( 1,2) ) 

GO TO 146 

145 IFtSin ,2) •NE,1. ) Al=TXl(l,IB+i)/Sin,2)/(51( I,2) + S2(l,2)) 

IF(S2( I ,2).NE.1.) A2=TX1(2, ie+l)/S2( I , 2 ) / ( S 1 ( I ,2 ) +S 2U , 2) ) 

146 CONTINUE 

A3=T( I, IB)/D1( I,2)/(D1(I,2)+D2{I,2)) 

A4=T(I, IB + 2)/02( I,2)/(Di( I,2)+D2( 1,2)) 

TNEW=( A1+A2+DXY«( A3+A4) » /A 
TNEW- t 1 ,-CME ( I , I B + 1 ) >DME=!=TNEW 
SUM-SUM+(T(I,IB+l)-TNEW)’l'(T( I , IB + 1 )-TNEW) 

T( I,I8+1)=TNEW 
150 CONTINUE 

C BOUND. POINTS LINES PARLL. TC Y AXIS .. LOWER AND TOP 
IF( ITER.EO.l ) WRITE(6,457) 

457 FORMAT (/,7X,« I « ,7 X , • Bl ' , 8X, • Cl • , 8X, * B2 * , 8X, * C2* , 8X , • B3* , 8X , • C3 * , 9X 
!,• A*,9X,»B*,9X,*C* ,/ ) 

DO 200 I=2,NX1 
XI=I-1 
B2 = XI’<'DX 
C LOWER 

C2=YUI) 

XM=(C2+0. OOOD/DY 

M=XM+2 

0M=X3(M) 

IF(I.GT.Nl) DM=X5(M) 

IF(ABS( B2-DM).LT.DX) GO TO 155 

B1=B2-DX 

Cl=Yl( I-l ) 

GO TO 156 

155 ai=DM ' 

C1=(FL0AT(M)-1.)>«'0Y 

156 CONTINUE 
DM=X4(M) 

IF(i.GT.NI) DM=X6(M) 

IF(ABS(D2-DM).LT.0X) GO To 160 

B3=B2+0X 

C3=Y1(I + M 
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160 

161 

165 

166 

170 

171 


174 

175 

176 

179 

180 

181 

490 


GO TO 161 

CONTINUE 

B3=0M 


C3 = FL0AT( M-1 )*0Y 

CALL SLC( B1 ,C1 ,B2 tC2,B3,C3 »A, B,C) 

IF( ITER. EG. 1) URITE{6,490)I»Bl,CltB2iC2»B3,C3rA,B,C 

XMN=B>2.’f'C’!=B2 

XT=FLOAT( ^'*-l )*DY 

IF(ABS(XHN).LE. 0.015) GO TC 170 

X=IC2-XT)’J=XMN-hB2 

IF(X.LE.B2) GO TO 165 

TTT = T( I,M}+(T( I + 1,M)-T( I,M) )'t(X-32)/DX 
GO TO 166 

TTT = T(I,M)-(T(I,M)-Tn-1,M) ) ( B2- X ) / DX 

XNL=(X-R2) >^(X-B2) + t Yl( I)-XT)«(Y1( I )-XT) 

XNL=SGRT( XNL) 

GO TO 171 
TTT=T(1tM) 

XNL=Xr-Yl ( I ) 


Al,-H( I , 1 ) /XK 

TNFW=( XNL*AL=i^TG( I , 1 ) +TTT ) / { 1 . +XNL L ) 

TNEW---I l.-OME)-=T( I , IN (1 )-l > + 0 ME’ 5 =TNEW 
SUM=SUM+(T( I,IN( 1 , I )-l )-TMEW )«(T( I, IN( 1 , I)-U-TNEW) 
TClf INU ,.n-l)=TMEVJ 
UPPER POINTS 
C 2 =Y 2 m 

XM=(C2+0.000l) /DY ' 

M=XM+l 

DM=X3(M) 

IFd.GE.N 2 ) n,M=X 5 (M) 

IF(ABS(B 2 -DM).LT. 0 . 0001 ) GO TO 174 
IFtABSI B 2 -DM).LT.DX) GO To 175 
B 1 =B 2 -DX 


ci=Y 2 d-n 

GO TO 176 
Bl = DM 

C1=FL0AT ( M-l )=«=DY 
0M=X4(M) 

IFd.GE.N2) DM=X6(M) 

IF(ABS(B2-DM).LT. 0.0001) GO TO 179 
IF { ABS( B2-DM) .LT.DX) GO TO 180 
R3=B2+DX 
C3=Y2( I+l ) 

GO TO 181 
B3=DM • 

C3 = FLQAT( M-1 >*DY • ' ' 

IFd.EQ.NXl) B3=X6 I IN{ 2, I )+l I 
IF(l.EO.NXl) C.3=FL0AT{ IN { 2d ) )=i=DY 
CALL SLC(Bl,CltB2.C2,R3fC3,A,B,C) 

TFdTER.EO.l) WPITE(6, 490) l,Bl»Cl,B2,C2,B3,C3,A,B,C 
F0RMAT(5X, 13 »2X,9F10.6 ) 

XMN=B+2.*C>>‘=B2 

IF(02d ,2 ) .EG. 1. ) GO TO 192 
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XT=FLOAT( M~l }-DY 

IF(ABS( XMN) .LE.O.OlSi GO TO 190 

X=(C2-XT)*XMN+D2 

IF(X.LE.B2 i GO TO 185 

TTT=T( I T{ I + l ,M)-T{ I » M ) ) ( X-B2 )/ DX 

GO TO 186 

185 TTT = T{ I ,M)-( T( I »M)-T( I -■ 1 » M ) ) « ( 32- X ) / DX 

186 XNL=( X-B2 )*(X-B2) + (Y2( I ) -XT ) - ( Y2 ( I ) -XT ) 

XNL=SQRT(XNL) 

GO TO 191 

192 X=DV*XMN+B2 

XNL = S0RT(nY’f=DY+(B2-X)>!^(B2-X) ) 

IF(X.LE-B2) GO TO 193 

TTT=T{ I ,M-1J+(T( 1+1,H-1) -T( I ,M-l) )*(X-B2)/DX 
GO TO 191 

193 TTT=T(I,M-1)<-IT{I,M-1)-T( I-I,M-1J )^!' ( B2-X) /DX 
GO TO 191 

190 TTT=T{I,.V) 

XNL=Y2(I)-XT 

191 AL=H(I,2)/XK 

TNEW=( XNL=>=AL’<'TG{ I,2)+TTT )/(l-+XNL*AL ) 
TNEW=U.-CKE)^T( I, IN(2 t I ) +1 ) +CME'-TNEW 
SUN=SUM* lTn»IN(2tI)+l )-TNEW )^(Tt I ,I M2 t I) + 1 )-TNEW1 
T( IfIN(2,I)+l)=TNEW 
200 CONTINUE 

C BOUND. POINTS.. NEARER AND FARTHER POINTS. 

C POINTS ON LINES PAPLL. TG X AXIS 
DO 300 J=2,NYi 
YJ=J-1 
B2=X3(J) 

C NEARER 

C2=YJ*0Y 

XM=(X3 { J)+0.000l) /nx 
N=XN+2 

IF{(Y2(N)-C2).LT.(C2-Y1{M) )) GO TO 220 

IF{ARS1C2-Y1(N) ).LT.ABS(C2-rL0AT(J-2)=«'DY) ) GO TO 206 
B1=X3{J-1J 
Cl=FLOAT{ J-2>*DY 
GO TO 208 

206 nl = FLOAT(N-l )=J'DX 
C1=Y1{N) 

208 IF(ABS( B2-FLnATlN-2 I'-'^nx) .LE. 0.0001 ) N=N-1 

IF(ABS{C2-FL0AT(J)^'CY).I.ToARS(C2~Yl(N-lJ )) GC TO 210 

B3=FL0ATI N-2 )-DX 

C3=Y1(Nt1) 

GO TO 240 
210 B3 = X3(JH) 

C3=FL0AT { J ) *DY 
GO TO 240 

220 IF(ABS{C2-Y2(N) ).l T.ARS(C2-FL0AT( J)*DYM gO TC225 
B1=X3(J>1) 

C1=FL0AT( J)=«--DY 
GO TO 226 
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225 ni = FLnAT(N-l)’<'DX 
Cl=Y2(N) 

226 IF(A0S(B2-FLaAT(N-2)’f'nX) .LT. 0-0001) N = M-L 
IF(A8.S(C2-PL0ATU-2)’?DY) .LT-ABS(C2-Y2{.Vl) ) ) GO TO 230 
B3 = FLnAT(N-2)’>’=0X 

C3=Y2(N-1) 

GO TO 240 
230 B3=X3(J-1) 

C3=FL0AT( J-2)*DY 
GO TO 240 

240 CALL SL0(Bl,CltB2,C2,B3,C3,A»n»C) 

IF ( ITER.EQ.l ) WRITS(6» 490) J,61,Cl,B2,C2, B3,C3,A, B,C 

XMN=B+2.*C'-'^02 

YT=FLOAT( N-1 )-''-DX 

IF(ABS( XMN) -GE.60, ) GO To 250 

Y=C2-( YT-B2) /XMN 

IF(Y.LE,C2) GO TO 245 

TTT=T( N,J ) + (T(N, J + 1 )-T {N, J) )=!=(Y-C2)/DY 
GO TO 246 

245 TTT=T(N, J)-(T(N, J)-T(M, J-1 ) )*tC2-Y)/DY 

246 XNL=(YT-B2)'^{YT-B2 ) + ( Y-C2 I Y-C2 ) 

XNL=5QRT{ XNL ) 

GO TO 251 

250 TTT=T(NtJ) 

XNL=YT-B2 

251 AL=HX( 1 t J)/XK 

TNEW=« TTT + AL=!'XNL=<'TGX( It J) ) / ( 1. »-XNL*AL) 

TNEW=I 1.-0ME>’-TX{ 1 , J) + GME*TNE'W 
SUM=SUM+(TNEW-TX( 1 ♦ J )}>>=( TNEW-TX ( 1 , J ) ) 

TX(1,J)=TNEW 
C FARTHER 

B2=X4(J) 

XN = { X4IJ) +o<.onoi ) /DX 
N-XNH 

!F (ABS( B2-FL0AT{ M-1)»!-DX) .LT. 0.0001) N=N-1 
IF((Y2(N)-C2).LT-(C2-Y1(N) ) ) GO TO 290 

IF( ABS( 02-FL0AT(N-1 )*DX) .LT.ABSt B2-X4U-1) ) ) GO TO 282 
BI--X41J-1 ) 

Cl=FL0ATU-2)-;'DY 
GO TO 283 

2B2. B1 = FL0AT( N-1 )--=DX 
CL=Y1(N) 

283 IF(ABS(B2--FL0AT(M)’S'DX).LT. 0,0001) N“N+l 

IFIABSI B2-FL0AT(M)*DX) ,LT*A05{B2-X4{ J-J-i) ) ) GG TO 284 
B3 = X4(j'-H) 

C3=FLGAT( J)-DY 
GO TO 261 

284 B3=FLOAT(N)*OX 
C3=YUN+1) 

GO TO 261 

290 IF ( AB5 ( B2-FL0AT (N-1 ) ^=DX) . LT. A8S ( D2~X4( J + 1) ) ) GO TO 292 
Bl=X4tJ^l) 

C1 = FL0AT( J )’!''DY 
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GO TO 293 

292 Bl=FLOAT(N^l )«DX 
Cl=Y2( N) 

293 IFiABS(82-FLOAT(N)’^DX),LT. 0.0001) N=N+L 

IF ( ABS( 82-FL0AT( N )=^DX) .LT. ABS ( B2-X4( J-l) ) ) GOTO 295 
B3 = X4( J-l ) 

C3 = FL0AT( J-2 j’f'DY 
GO TO 261 

295 B3 = FL0AT(M)=^DX 
C3=Y2(N+1) 

261 CALL SLOIBl ,CltC2,C2,B3,C3,A,a,C) 

IF( ItER.EO.I ) WRITE(6,49Q)J,Dl»Cl,B2,C2»B3iC3tA,B,C 
XMN=B + 2.=:«C«fl2 

, • YT = FL0ATI N-1 )*DX 

IF(ABS(XMN).GE.60. J GO TO 270 
Y=C2-( YT-B2) /XMN 
IF(Y.LE.C2J go to 265 

TTT=T( N,J ) <-(T( N, J + 1 )-T (N, J ) )=«'(Y-C2)/DY 
GO TO 266 

265 T r T = T ( N , J ) - { T { N » J ) -T ( N t J- 1) ) "M C2- Y ) / DY 

266 XNL = {YT-B2)»^(YT-B2 )*-(Y-C2)^MY-C2) 

XNL=SORT( X\'L) 

GO TO 271 

270 TTT=T(NtJ) 

XNL=R2~YT 

271 AL=HX(2,U )/XK 

TNEW = ( TTT + AL'-:'XNL=:=TGX (2 , J I ) /( l.+XML^AL) 

TNEW={ l.-0r^E)=!=TX(2, J) + CME*TNEW 
SUM=SUM+(TNEW-TX(2, J) ) ’!= ( TNEW-TX { 2, J ) ) 

300 TX(2,J)=TNEW 

C BOUNDARY POINTS FOP J=N4 TG N5 

IF(N4.EQ.99) GO TO 351 
DO 350 J=N4,N5 

YJ=J~1 

B2=X5(J) 

C2=YJ*DY 

XN=(X5(J) +0.0001 )/DX 
N=XN+2 

B1 = FL0AT( N-1 J+'OX 
C1=Y1(N) 

B3=B1-DX 

C3=Yl(N-X) 

IF{B3.EQ.B2> B3--=Q3-0X 
IF(C3,E0.C2) C3=YUN-2) 

CALL SL.QT Bl,CI,B2,C2tB3»C3f A,B,C) 

XMN=B+2.vC^B2 

IFf ABS{XMN).Gr-.60. ) GO TO 320 

YT = FLOAn J )'-DY 

X--=H2+(G2-YT) 

IND=X/DX 

INl=IND+l 

TTT=T( IND . J+l) + ( T n Nl, J-+1 )-T ( I NO, J + 1 ) ) >!' ( X- F LOAT { IN D- 1 ) * DX ) / OX 
XNL={X-B2) X-B2) + ( VT-C2)’^-(YT-C2 ) 
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XNL=SQRT( XNL) 

GO TO 322 
320 TTT=T(N,J) 

XML=FL0A-nM-i)>!‘DX-X5U ) 

322 AL=HXi( 1 t J ) /XK 

TNEW=( TTT+AL^-XNL=:'TGXll 1, J) ) / ( I« f AL*XNL ) 
TNEW=U .-GME)*TXL( I , J ) +QM.5«TNEvJ 
SUM=SUM+{ TXl(l»J)“TNEW)=!'-(TXi ll,J)-TNEW) 

TXK L, J)=TNEW 
B2=X6IJ) 

XN=(X6( J)-5-0.000l)/DX 
N=XN+1 

IF (FLOAT(N)^OX.EQ„D2) N=M-1 

IF {ABS{Y1(N)-C2) .LT.DY) GO TO 325 

IF(ABS( FLOAT(M)-'.''OX-82).LT.O.OOOU GO TO 324 

IF(N.GE.NX) GO TO 327 

B1 = FL0AT{N)^.= DX 

Cl=Y2(N+l) 

GO TO 323 
327 B3=X6(J-l) 

C3=FLGAT{ J-2 )=^-DY 
GO TO 323 

324 B1=FL0AT( N+1 )-0X . 

C1=Y2( M+2 ) 

GO TO 323 

325 Bl=FLOAT{N-l)*DX 
CI.=Y1{N) 

323 IF(ABSCY2(N)-C2).LT.DY) GO TO 326 
B3-=X61J+I) 

IF ( ( J + n*GT.N5) B3 = X4(J+1) 

C3 = rL0ATt J )^'DY 
GO TO 328 

326 B3=FL0AT(N-1)=^=DX 
C3=Y2(N) 

323 CONTINUE 

CALL SLO< BltCl ,B27C2,B3,C3 tA,B,C) 
XMN^B+2.>:=C^=B2 
IFU.EQ.N6) GO TO 335 
XT = B3 

Y=C2-(XT-B2)/XMN 

TTT = T(N-l » J)-{ T(N-1 , J)-T(N-l » J-L) } ’!'( C2-Y1 /OY 
XNL={C2-J)*(C2-J)<-IB2^B3)*(B2-B3) 

XNL=SQRT( XNL) 

GO TO 336 

335 TTT^TI N.XL ,M6) 

XNL=B2-FL0AT (NX-2)*CX 

336 AL=HX). {2, J)/XK 

TMEW=(TTT + AL'i^XNL*TGXl( 2, J))/ (l. + AL=i=XNU 
SUM--S'JM+(TX112, J)-TNFW)v( vxi (2,J )-TNEW) 
TNHW=( l.-nME)*TXL( 2, J) -t-OYE*TNEW 
TXI(2, J) = TilPW 

350 CONTINUE 

351 CONTINUE 


35 
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IV G LEVEL 21 MAIN DATE = .74273 15/22/49 

T( 1,4»=TX( 1,4» 

TC1,6)=TX(1,6) 

T(40t4)=1982« 

T(40,5)=TX1(2,5) 

TC40,6)=i982. 

WRITE(6,420) ITER, SUN 

420 FORMAT (4X, • ITERATION = • , I 3 ,4X , * SUK= • , E12 .5 ) 

IF(SUM.LE.SUMM) GO TO 400 
IFMTER/NTE^NTE.BQ, ITER) GOTO 400 
SUM=0. 

ITER=ITER+1 
GO TO 125 
400 CONTINUE 

WRITE(6,458) 

IF(SUM.LE.SUMM) WRITFC6,46i) 

461 FORMAT!/, 20X,' THE FINAL ITERATIVE* ) 

440 F0RMAT(3X,I2,i2F10.3/) 

WRITE(6,455) 

DO 450 1=1, NX 

WRITE(6,440) I , ( T( I , J) , J = l ,NY ) 

450 CONTINUE 

IF! SUM.LE.SUMM) GO TO 410 
IF! ITER.GE.NITE) GO TO 410 
415 CONTINUE 

ITER=ITER+1 

SUM=0, 

GO TO 125 
410 CONTINUE 

WRITE! 6,460) 

460 FORMAT! IHl ,33X, • ISOTHERMAL LINE LOCATIONS » / ,9X, • I • ,4X, * J* ,3X, *T • , 
14X, 'T-LOW* ,6X,*T-HIGH* , 8X , * F P AC* , / ) 

DO 500 1=2, NXl 

IA=IN!1,I)-1 

IB=IN!2,I) 

IT1=T! I ,IA ) 

IT2=T! I , 18+1 ) 

IF! IT1,GE.IT2) GO TO 470 

IL=IT1 

IH=IT2 

GO TO 471 

470 IL=IT2 
IH=IT1 

471 CONTINUE 

on 482 II=TL, IH 
DO 480 J=I A, IB 

IFmi,J).LT.FLnAT!IT).AND.T!l,J + l).GT.FLnAT!im GO TO 48 1 
IF!T! I, J) .GT.FLOAT! in.AND.T! I,J + i).LT.FLOAT! II) ) GO TO 48 3 
GO TO 480 

481 RX=!FLaAT!II)-T(I,J))/!T( I,J+l>-T!I,J)) 

GO TO 484 

483 RX=!T! I ,J )-FLOAT! 1 I ) )/!T( I ,J)-T! I , J+in 

484 CONTINUE 

WRITE!6*485) I , J , 1 1 , T ! I , J ) ,T ! I , J + 1 ) , RX 


36 



iv G LEVEL 2 1 


MAIN 


DATE = 74273 


15/22/49 


• w 


485 F0RMAT(5X,3I5,3Fll.4> 

480 CONTINUE 
482 CONTINUE 
500 CONTINUE 

CALL BUN(T»NX,NY) 

505 CONTINUE 
STOP 
ENU 


U 
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IV G LEVEL 21 SLO DATE = 74273 

SUBROUTINE SLR ( XL , Y I , X2t Y2 ,X3 , Y3 t A, B , C ) 

C FINDS EQUATION THRU THREE f»OINTS AS Y= fi + R$X+C=«'X*X 
CALL UNDFLW 

IF(A8S( Xl-X2).LE.O.OOl.ANU.A8S(X2-X3).LE.O.OOl) GO TO 10 

XS3=X3*X3 

XS2=X2=!‘X2 

XSl=Xl«Xl 

Dl=X2*XS3”X3'<'XS2 

D2 = XI^XS3-X3<'XS1 

D3=X1*XS2-X2'‘XS1 

D=D1-D2+D3 

IFCABS(D).LE. 0.000001) GO TO 10 
A1 = Y1*( X2^ XS3-X3’«=XS2) 

A2=Y2’<‘( Xl>i=XS3-X3’i'XSl) 

A3=Y3=^( Xl*XS2-X2^XSl) 

A=(A1-A2+A3)/D 
B1“Y2«XS3“Y3*XS2 
B2=Yl«XS3-Y3vXSl 
B3sYl«XS2-Y2’i'-XSl 
B=(BI-B2+B3)/0 
C1 = X2’!^Y3-X3*Y2 
C2 = Xl’t'Y3-X3*Yl 
C3=X1«Y2~X2^'Y1 
C=(C1-C2+C3)/D 
RETURN 
10 A=100. 

B=100. 

C=100. 

RETURN 

END 


15/22/49 
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IV G LEVEL 21 BUN DATE “ T4273 

SUBROUTINE 8UNIT,NX*NY) 

DIMENSION T(40,20l 

LOGICAL*! KP(60 j,LlT(l5)/»0* , * I* * *2* ,*3» , f *5* ♦ *6* , 
**=*/ 

TMA-Td ,1) 

TMI-T{1,1) 

DO 10 1=1, NX 
DO 10 J=1,NY 

. IF1T( I,J).GT.TMA) TMA=Tn,J) 

IFITU, J) .LT.TMI ) TKI = T(I,J) 

10 CONTINUE 

WRITE(6,20J TMA,TMI 

20 . F0RMAT(1HI,5X,*MAX=*,F12.5,2X,'M1N=* , FI2 . 5/ , 5X , * TE M P. 

1 : 0-l-2-3-4-5-6-7-8-9“.-?-+-$-*» , 5(/)I 
WRITE(6,5) 

5 F0RMAT(35X, ‘PICTORIAL VIEW OF TEMPERATUR E » , 15( / I ) 

RC = 14.99«599/(TMA-TMI ) 

DO 30 JX=I,NY 
J=NY+1-JX 
DO 29 1=1, NX 
S=RC*(T( I, JI-TMI ) 

29 KPU )=LIT(l + IFIXlSn 
WRITE! 6,35) (KP ( I ) , 1=1 ,NX ) 

30 CONTINUE 

35 FORMATdH , • . ' , 25X , 50( A 1 , 1 X ) ) 

RETURN 

END 


ORIGlNAEPAlSBi^ 
OF POOR QUALTH] 


15/22 /H9 

7» , ‘S* , *9* ,« 
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APPENDIX C 


SAMPLE PROGRAM OUTPUT 

The program output includes the boundary point coordinates 
which can be used to check the program input parameters. The 
program also prints out, e, for every iteration and this allows 
monitoring of the convergence process. After the final 
iteration, the grid temperatures are printed as well as the 
isothermal lines locations. A sample of the program output 
is given in the pages 41 through 56 . 
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TEMPERATURE DISTRIRUTION IN A TURBINE BLADE 


DX - C„0500 OY = 0.0250 OMEGA = 1.600 THERM. COND. = 1.0000 

NX = 40 NY = 12 MAX. ERR.= 0.0200 




IM.’tC;'! IS TMf PlUNT HN TMf rtrONI'APY AN!5 I (• 1 . C I I , ( 'ii t C3 I A« t IHt SUPPDUNOINC. I’UINTS UN hlADf 
A.tiiC ARf THl CMIFFICIFUTS OF DIE I’AFAElllCC. CUI'Vf. Y« AtB* X«C • **X 


I 


Cl 

82 

C2 

83 

C3 

A 

B C 

2 

0. 04 0000 

0.050000 

0.050000 

0. 045000 

0. 100000 

0.028000 

0.075333 

-0.73999B 2.666591 

2 

0.028000 

0. 15U000 

0.050000 

0. 163POO 

0. 100000 

0. 1 801)00 

0.128516 

0.862729 -3.484890 


0, 09 0000 

0.045000 

0 . lonotn 

0. 028000 

0. 150000 

0.021000 

0.077000 

-0.640001 1.99999B 

3 

0. 08 0000 

0, 17'joon 

0. 100000 

0. 1 801)00 

0. 1 50000 

0. 193000 

0.156143 

0.2242P6 0.142830 

A 

0.120000 

0. 02 5000 

0. 150|)00 

0. 021000 

0. 200000 

0. 025000 

O.OUU99U 

-0.053307 2.66657B 

4 

0. 100000 

0. 180000 

0.150000 

0. 193000 

J. 200000 

0.203000 

0.146001 

0.409980 -0.599907 

S 

0. 150000 

0.021000 

0.200000 

0.02 5001 

0.2 500 00 

0. 035000 

0.045000 

-0.339995 1.199900 

5 

0. 185000 

0. 200000 

0. 200000 

0.203000 

0. 2 50000 

0. 21 1000 

0. 140737 

0.636956 -0. 615390 

& 

0.200000 

0.025000 

0. 250000 

0.035000 

0.300000 

0.045000 

-0.014999 

0.199996 0.000015 

A 

0. 200000 

0.203000 

0.250000 

0.21 1000 

0.3000 00 

0.220000 

0.180999 

0,069991 0.200048 

7 

0. 25 0000 

0.035000 

0. 300000 

0,045000 

0. 32 0000 

0.050000 

0.03H572 

-0.192072 0.714225 

7 

0.250000 

0.211000 

0. KlUOOO 

0.220000 

0.350000 

0.22 7000 

0.136001 

0.399985 -0.400119 

B 

0.300000 

0.045000 

0.350000 

0. 055000 

0.400000 

0.068000 

0.047999 

-0.190004 0. 599594' 

B 

0. 340000 

0.225000 

0. 350000 

0.227000 

0.400000 

0.233000 

-0.001676 

1,120079 -1.333167 

9 

0.350000 

0.055000 

0. 400000 

0.06 3000 

0.440000 

0.076000 

-O.U.8126 

0.967934 -0.944040 

■5 

0. 350000 

0.227000. 

0. 400000 

0.233000 

0.4 50000 

0. 238000 

0.156997 

0.270025 -0.200045 

10 

0.400000 

0. 068000 

0.450000 

0.077000 

0.500000 

0.086000 

-0.004005 

0.180019 0.0 

10 

0.900000 

0.233000 

0.450000 

0.238000 

0.500000 

0,243000 

0. 193003 

0.100009 0.000238 

11 

0.450000 

0. 077000 

0. 500000 

0. 086000 

0. 55000U 

0.095000 

-0.003984 

0.179961 0.0 

11 

0. 950000 

0.23BOOO 

0. 500000 

0.243000 

0.550000 

0.245000 

0.058036 

0.670009 -0. 600000 

12 

0.500000 

0. 086000 

0.550000 

0.095000 

0.575000 

0. lOOOOO 

0.069342 

-Q. 099905 0.266479 

12 

0. 500000 

0.243000 

0. 550000 

0.245000 

0.600000 

0.249000 

0.333022 

-0.380336 0.400429 

13 

0. 56 0000 

0.095000 

0.600000 

0. 105000 

0.650000 

0. 113000 

-0. 146953 

0.659878 -0.400060 

13 

0. 550000 

0.245000 

0.600000 

0.249000 

0.650000 

0.251000 

0.069040 

0.539526 -0.400000 

lA 

O.oOOOOO 

0. 105000 

0.650000 

0. U3000 

0. 700000 

0.120000 

-0.069015 

0.410035 -0.200048 

lA 

0. 630000 

O.2SO00Q 

0,650000 

0.251000 

0.700000 

0.254000 

0.276955 

-0.131915 0.142979 

13 

0.650000 

0. 113000 

0. 700000 

0. 120000 

0.735000 

0.125000 

0.037314 

0.094226 0.034054 

13 

0.650000 

0.251000 

0. 700000 

0.254000 

0. 750000 

0.255000 

0.030036 

0.600000 -0.400000 

IF 

0. 700000 

0. 120000 

O.75U000 

0. 128000 

0.800000 

0.133000 

-0. 307027 

1.030043 -0.600381 

16 

0.700000 

0.254000 

0. 750000 

0.255000 

0.800000 

0.255000 

0. 134984 

0.309967 -0.199809 

17 

0. 750000 

0. 128000 

0. 800000 

0. 133000 

0. 850000 

0. 138000 

0.052992 

0.099905 -0.000238 

17 

0. 75 0000 

0.255000 

0.800000 

0.255000 

0.850000 

0.255000 

0.254992 

0.0 0.0 

IB 

0. BOOOOO 

0. 133000 

0. 850000 

0.138000 

0.900000 

0.143000 

0.052977 

0.099905 -0.000238 

IB 

0.800000 

0.255000 

0.850000 

0.255000 

0.900000 

0.254000 

0. 118979 

0,329995 -0.200286 

19 

0.850000 

0. 13B00Q 

0. 900000 

0. 143000 

0.950000 

0.145000 

-0.406026 

1.150214 -0.600143 

19 

0. 850000 

0.255000 

0.900000 

0,254000 

0.950000 

0.25100U 

-0.034022 

0.680257 -0.400095 

20 

0. 900000 

0.143000 

0. 950000 

0. 14SJ00 

1. 000000 

0. '48000 

0.277950 

-0.329917 0.200000 

20 

0. 900000 

0.2 54000 

0.950000 

0.251000 

1.000000 

0.250000 

0.649940 

-0.799762 0.399762 

21 

0. S50000 

0.145000 

1. 000000 

0. I4S000 

1.040000 

0.150000 

-0.016815 

0.275305 -0.110452 

21 

0. 950000 

0.25100C 

1.000000 

0. 250000 

1.049999 

0.248000 

0.080087 

0.370026 -0.200143 

22 

0.999999 

0. 148000 

1. 049999 

0. 161000 

1.095998 

0.154000 

0.087965 

0.059876 -0.000239 

22 

0.999999 

0.250000 

1.049999 

0.248000 

1.099998 

0,245000 

0.079914 

0.370229 -0.200143 

23 

1.099999 

0.15100Q . 

1. 099999 

0. 154000 

1.149999 

0.156000 

-0.372297 

0.916825 -0.398764 

23 

1.099999 

0.248000 

1.099999 

0.245000 

1.149999 

0,242000 

0.310792 

-0.059886 -0.000238 

2A 

1.099999 

0.154000 

1.150000 

0. 155000 

1. 199999 

0.155000 

-0.120099 

0.468631 -0.199144 

2A 

1.099999 

0.245000 

1.150000 

0. 242000 

1.199999 

0.238000 

0.058564 

0.388783 -0.199144 

2fi 

1. 150000 

0, 1.55000 

1. 200000 

0. 155000 

1.2499 99 

0.155000 

0.154988 

0.0 0.0 

25 

1.1 50000 

0.242000 

1.200000 

0.238000 

1,249999 

0.232000 

-0.216537 

0.857414 -0.398527 

2(i 

1.199999 

0.155000 

1.249999 

0. 155000 

1.299998 

0.154000 

-0.145172 

0.490219 -0.200143 

2fc 

1.199999 

0. 238000 

1.249999 

0.232000 

1.299998 

0.228000 

0.982362 

-1.100429 0.400286 

27 

1.299999 

0. 155000 

1.299999 

0. 154000 

1.349993 

0. ISIOOO 

-0.468007 

0.996911 -0.3'J3764 

27 

1.299999 

3.232000 

1.299999 

0.228000 

:.,3.'3ooa 

0.225000 

-0.5293 17 

1.271464 -0.529040 

2B 

1.299999 

0. 15400Q 

1. 349999 

0. ISIOOO 

•1.399999 

0.147000 

-0,118138 

0.468631 -0.199620 

28 

1.299999 

0.228000 

1.349999 

0. 222000 

1, 399999 

0.215000 

0.033582 

0.408745 -0.199382 

29 

1,399999 

0.151000 

1.400000 

0.147000 

1.449999 

0.142000 

-0.119170 

0.470420 -0.200143 

29 

1.399999 

0.222000 

1.400000 

0.215000 

1.449999 

0.208000 

0.411051 

-0.140267 0.0 

30 

1,900000 

0.147000 

1.450000 

0. 142000 

1. 499999 

0.138000 

0.691154 

-0.667538 0.199144 

30 

1.900000 

0.215000 

1.450000 

0. 208000 

1.499999 

0.200000 

0.005629 

0.420945 -0.199144 

31 

1.949999 

0.142000 

1.499999 

0. 138000 

1. 54999B 

0.133000 

-0.177242 

0.510496 -0.199905 

31 

1.449999 

O.20S00O 

1.499999 

0.200000 

1.549998 

0.193000 

0.875194 

-0. 750000 0.200143 

32 

1.499999 

0. 133000 

1. 549999 

0. 133000 

1.599998 

0. 127000 

-0.176034 

0.508317 -0.199144 

32 

1.499999 

0.200000 

1. 549999 

0.193000 

1- 599998 

0.183000 

-0.901226 

1.684173 -0.598384 

33 

1.549999 

0.133000 

1.599999 

0. 127000 

1.649999 

0.123000 

1.311604 

-1.360725 0.400048 

33 

1.549999 

0. 193000 

1. 599999 

C. 183000 

1.645000 

0. 175000 

1.083652 

-0.937500 0.234096 

34 

1.625000 

0.125000 

1.650000 

0. 123000 

1.699999 

0. 1 1 5000 

-2.585740 

3.389520 -1.058712 

34 

1,599999 

0. IS’OOO 

1.650000 

0, 174000 

1.699999 

0,165000 

0.471091 

-0.180105 -0.000239 

35 

1 .650000 

0. 123000 

1. 7000 00 

0.115000 

1.749999 

0.108000 

0.948414 

-0.830153 0.199905 

35 

1.650000 

0.174000 

1.700000 

0.1650G0 

1.749999 

0,153000 

-1.212815 

1.831106 -0.600429 

36 

1. 699999 

0.115000 

1. 749999 

0. 108000 

1. 799993 

0. lOOOOO 

-0.242292 

0.550334 -0.200143 

36 

1.699999 

0.165000 

1.749999 

0. 153000 

1, 7 70000 

0,150000 

4.363598 

-4.637669 1.275337 

37 

1.749999 

0. lOSOOO 

1.799999 

0.100000 

1,849998 

0.097000 

3.526705 

-3.698194 0.996673 

37 

1.749999 

0.153000 

1. 799999 

0. 143000 

1.849998 

0,133000 

0.501708 

-0.199144 -0.000238 

3B 

1.799999 

0. t 00000 

1.849999 

0.097000 

1.899999 

0.085000 

-5.789 301 

6.513358 -1.801049 

38 

1,799999 

0, 143000 

1.849999 

0. 133000 

1.8 84999 

0.126000 

-0.616411 

1.026842 -0.335737 

39 

1.049999 

0.097000 

1.900000 

0. 085000 

1.940000 

0.100000 

24.526886- 

■25.830017 6.824074 

39 

1.849999 

0.133000 

1 .900000 

0. 120000 

1.884999 

0.125000 

-6.883969 

7.743055 -2.134259 

2 

0.150000 

0, 021000 

0. 120000 

0. 02 5000 

0, 100000 

O.O2BO0O 

0.047001 

-0.723317 0.333261 

2 

0.150000 

0. 021000 

O.2UOUO0 

0.02 5000 

0.2 50000 

0.035000 

0.045000 

-0.339995 1.199980 

3 

O.OOOOOO 

0. 045000 

0.040000 

0.050000 

0.010000 

0.075000 

0.086667 

-1.249997 8.333333 

3 

0.300000 

0.045000 

0. 320000 

0. 050000 

0. 350000 

0.055000 

-0.190000 

1.283372 -1.666832 

4 

0.040000 

0.050000 

O.OIOOOO 

0.075000 

0.0 

0. lOOOOO 

0. lOOOOO 

-2.916664 41.566534 

4 

0.400000 

0.068000 

0.440000 

0.075000 

0.450000 

0.077000 

0.085938 

-0.243860 0.497774 

5 

O.OIOOOO 

0.075000 

0.0 

0 . lonooo 

0.005000 

0.125000100.000000100. 00 0000 1 00 . 000000 

5 

0. 550000 

0. 095000 

0.575000 

0. lOUOOO 

0.600000 

0.105000 

-0.014836 

0.199643 0.000119 

6 

0. 02 8000 

0. 150000 

0.005000 

0.125000 

0.0 

0 . 1 ooooo 

0. lOOOOO 

5. 698756 *••••••**• 

C 

0.70QOOO 

0.120000 

0. 735000 

0. 12 5000 

0. 750000 

0.128000 

0.606264 

-1.492646 1.13800B 

7 

0.050000 

0.163000 

0.028000 

0. 150000 

0.005000 

0.125000 

0.118022 

1.450722-11.023401 

7 

l.OOQOOO 

0.148000 

1 . 040000 

0. 150000 

1.049999 

0.151000 

1.136719 

-1.988094 1.000000 

8 

0. lOOOOO 

0.180000 

U. 080000 

0.176000 

0.050000 

0.163000 

0.130999 

0.790001 -3.000001 

8 

1. 599999 

0.18 3000 

1.645000 

0.175000 

1.650000 

0. 174000 

-0.766335 

1.357954 -0.471591 

9 

0.200000 

0.203000 

.0. 18 6000 

0.200000 

0, 150000 

0. 193UU0 

0.1 63001 

0.200255 -0.000142 

9 

1.450000 

0.208000 

1.500000 

0.200000 

1.549999 

0.193000 

0.877814 

-0 . 753113 0.2 0 0 6 70 

10 

0. 350000 

0.227000 

0. 340000 

0.22 5000 

0.300000 

0.270000 

0.335491 

-0.834979 1.498975 

10 

1.299999 

0.228000 

1.325000 

0. 225000 

1. 349999 

0.227000 

0.382813 

-0.117424 0.001894 

11 

0,650000 

0.251000 

0.630000 

0. 250000 

0.600000 

0.249000 

0.355166 

-0.377734 0.333996 

11 

0. 950000 

0.251000 

l.OOOUUU 

0.250000 

1.049999 

0.248000 

0.080157 

0.369700 -0.199857 





ITERATION 

= 

L 

SUM^ 0.L7129E 

05 

ITERATION 


2 

SUM= 0.10122E 

05 

ITERATION 

- 

3 

SUH= 0.477B66 

04 

I TERATION 

— 

4 

SUM= 0.23836E 

04 

ITERATION 

= 

5 

SUM= 0.205005 

04 

ITERATION 


6 

SUM= 0.16037E 

04 

ITERATION 

= 

7 

SUM= 0.14074E 

04 

ITERATION 


a 

SUM= 0- 115165 

04 

ITERATION 

= 

9 

5UM= 0.11^135 

04 

ITERATION 

= 

LO 

SUM= 0.96032E 

03 

ITERATION 

= 

11 

SUM= 0.91320E 

03 

ITERATION 

= 

12 

SUN= 0.81725E 

03 

ITERATION 


13 

SUM= 0.75635E 

03 

ITERATION 

= 

14 

SUM= 0.69111E 

03 

ITERATION 

= 

15 

Sl)iM= 0.64236E 

03 

ITERATION 

=: 

16 

SUM= 0.587745 

03 

ITERATION 


17 

SUM= 0.54829E 

03 

ITERATION 

- 

10 

SUM= 0.50517E 

03 

ITERATION 


19 

SUM- 0o46995E 

03 

ITERATION 

=: 

20 

SUM= 0.43643E 

03 

ITERATION 


21 

SUM= 0.40629E 

03 

ITERATION 

= 

22 

SUM= 0.37849E 

03 

ITERATION 

= 

23 

SUM= 0.353625 

03 

ITERATION 


24 

SUM= 0.32970E 

OS- 

Iteration 

= 

25 

SUM= 0.30914E 

OS 

ITERATION 

■= 

26 

5UM= 0.2892 IE 

03 

iteration 


27 

SUM= 0.27139E 

03 

ITERATION 

= 

28 

SUM= 0.25468E 

03 

iteration 


29 

SUM= 0.240005 

03 

ITERATION 


30 

SUM= 0.22515E 

03 

Iteration- 

= 

31 

SUM= 0.21254E 

03 

iteration 

= 

32 

SUM= 0. 19997E 

03 

ITERATION 

= 

3 3 

SUM= 0.13928E 

03 

iteration 


34 

SUM= 0.17B43E 

03 

ITERATION 


35 

SUM= 0.168795 

03 

iteratigm 


36 

SUH= 0.15971E 

03 

ITERATION 

= 

•37 

SUM= 0. 15140E 

03 

HERAT ION 

= 

38 

SUM= 0.14321E 

03 

iteration 

r: 

39 

SUM= 0. 13610E 

03 

ITERATION 

= 

40 

SUM= 0.1288SE 

03 

HERAT ION 

— 

41 

SUM= 0. 12272E 

03 

ITERATION 

= 

42 

SUM= 0.116445 

03 

ITERATION 

= 

43 

SUM= 0.110795 

03 

ITERATION 


44 

SU.M= 0.10 526 5 

0 3 

ITERATION 

= 

45 

SUM= 0.100315 

03 

ITERATIGM 

= 

46 

SUM= 0.95471E 

02 

ITERATION 

= 

47 

5UM= 0.90099E 

02 

ITERATION 

= 

48 

SUM= 0. 86 74 IE 

02 

ITERATION 

= 

49 

SU.M= 0.8265 IE 

02 

ITERATION 

= 

5 0 

SUM^ 0.73813B 

02 
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TEMPERATURE AT GRID POINTS 


1 

1 

2 

3 

4 

5 

6 

7 

a 

9 

10 

11 

12 

1 

1961.000 

1961.000 

1961.000 

1953.493 

1952.095 

1954.163 

1961.000 

1961.000 

1961.000 

1961.000 

1961.000 

1961.000 

2 

1966.000 

1952.318 

1951.977 

1962.025 

1952.434 

1953.486 

1955.266 

1956.594 

1998.500 

1998.500 

1998.500 

1998.500 

3 

1975.500 

1952.739 

1952.077 

1952.193 

1952.566 

1953.208 

1954.050 

1954.897 

1955.272 

1998.500 

1998.600 

1998.500 

A 

1952.426 

1952,111 

1952.009 

1952.098 

1952. 336 

1952.874 

1953.543 

1954.387 

1955.243 

1993.350 

1998.350 

1998.350 

5 

1975.200 

195 2.773 

1951.976 

1951.994 

1952. 193 

1952.575 

1953.105 

1953-743 

1954,423 

1954.785 

1998.200 

1998.200 

6 

1976.200 

1952. 551 

1951.949 

1951.958 

1952. 125 

1952.468 

1952.997 

1953.723 

1954,665 

1955.395 

1998,150 

1998.150 

7 

1977.200 

1952.642 

1952.284 

1952.112 

1952. 157 

1952. 414 

1952.085 

1953,593 

1954,582 

1955.73 7 

1998. 100 

1996. 100 

a 

1977,900 

1977. 900 

1952.967 

195 2.22 3 

1952.177 

1552.309 

195 2. 622 

1953.097 

1953.693 

195 4. 304 

1954.637 

1997.500 

9 

1973.600 

1978.600 

1952.995 

1952.574 

1952.313 

1952.306 

1952.518 

1952.933 

1953.536 

1954.320 

1954.890 

1996.900 

to 

1973.600 

1978.800 

1978,800 

1953.180 

1952,294 

1952.243 

1952.423 

1952.819 

1953.433 

1954,263 

1954.990 

1996.700 

11 

1979.000 

1979.000 

1979.000 

1952,812 

1952. 145 

1952.106 

1952.271 

1952.646 

1953.234 

1954.039 

1954.877 

1996.500 

12 

1979.000 

1979.000 

1979.000 

1952.293 

1952,050 

1951.992 

1952.105 

1952.419 

1952.948 

1953.706 

1954.576 

1996.500 

13 

1979.000 

1979.000 

1979.000 

1979.000 

1952.889 

1952.03 7 

1951.971 

1952-136 

1952.532 

1953.184 

1954.134 

1996.500 

!A 

1973.900 

1978.900 

1978.900 

1978.900 

1952.518 

1951.881 

1951.754 

1951.762 

1951.905 

1952.154 

1952.434 

1952.719 

15 

1978.800 

1978.800 

1978.800 

1978.800 

1952.616 

1952.211 

1951.858 

1951.700 

1951.705 

1951.843 

1952.070 

1952.472 

16 

1978.400 

1978.400 

1978.400 

1978,400 

1978.400 

1953.045 

1952.056 

1951.847 

1951.840 

1952.027 

1952.410 

1952.858 

17 

1978.000 

1978.000 

1978.000 

1978.000 

1978.000 

1953. 102 

1952.256 

1952,057 

1952.065 

1952.276 

1952.701 

1953.158 

IS 

1976.900 

1976.900 

ig76c«00 

1976.900 

J576.900 

1953.136 

1952.462 

1952.313 

1952.324 

1952.507 

1952.911 

1953. 372 

19 

1975. 800 

1975.800 

1975.800 

1975.800 

1975.800 

1953.324 

1952.832 

1952.803 

1952.748 

1952.700 

1952.781 

1953.206 

20 

1974.850 

1974.850 

1974.850 

1974.850' 

1974.850 

1954.561 

1954. 164 

1954.240 

1953.993 

,1953.241 

1951.490 

1951.804 

21 

1973.900 

1973.900 

1973.900 

1973.900 

1973.900 

1957.775 

1957-516 

1957.746 

1958.091 

1958,703 

1959.905 

1993.500 

22 

1973.050 

1973.050 

1973.050 

1973.050 

1973.050 

1973.050 

1961.964 

1961.430 

1961.720 

1962.350 

1963.259 

1992.500 

23 

1972.200 

1972.200 

1972.200 

1972.200 

1972.200 

1972.200 

.1964.589 

1964.200 

1964.502 

1965.118 

1965.892 

1991.500 

24 

1971.700 

1971.700 

1971.700 

1971.700 

1971.700 

1971.700 

1966.511 

1966.190 

1966.494 

1967.083 

196 7.724 

1990,050 

25 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1967.716 

1967.465 

1967.779 

1968.332 

1968.830 

1988.600 

26 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1968. 44B 

1968.229 

1968.516 

1968,977 

1969.291 

1986.900 

27 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1968.870 

1968.667 

1968.955 

1969.393 

1969.587 

1985.200 

2fi 

1970. 900 

1970. 900 

1970.900 

1970.900 

1970.900 

1970.900 

1968,947 

1968.775 

1969.177 

1969.767 

1983.350 

1983.350 

29 

1970.600 

1970. 600 

1970.600 

1970.600 

1970.600 

1967.701 

1967.559 

1968. 160 

196 6.719 

1969.154 

1981.500 

1981.500 

30 

1969.650 

1969.650 

1969.650 

1969.650 

1969.650 

1967.812 

1967.677 

1968.029 

1968.427 

1968.678 

1979.550 

1979.550 

31 

1968. 700 

196 8.700 

1968.700 

1966.700 

1968.700 

1967.801 

1967.692 

1968.010 

1968.461 

1977.600 

1977.600 

1977.600 

32 

1968.350 

1968.350 

1968.350 

1968.350 

1968.350 

1967.504 

1967.411 ■ 

1967.693 

1968.015 

1975.800 

1975.800 

1975.800 

33 

1968.000 

1968.000 

1968.000 

1968.000 

1968.000 

1966.953 

1966.863 

1967.124 

1967.305 

1974.000 

1974,000 

1974.000 

34 

1967.250 

1967.250 

1967.250 

1967.250 

1965.681 

1965.583 

1966.057 

1966.424 

1972.500 

1972.500 

1972.500 

1972.500 

35 

1966.500 

1966.500 

1966.500 

1966.500 

1965. 67H 

1965.598 

1965.867 

196 6.046 

1970.500 

1970.500 

1970.500 

1970.500 

36 

1965.750 

1965.75 0 

1965.750 

1965.750 

1965.290 

1965.337 

1966.127 

1966.188 

1969.000 

1969.000 

1969.000 

1969.000 

37 

1965.000 

1965.000 

1965.000 

1965.000 

1963.041 

1962.869 

1963.242 

1967.500 

196 7.500 

1967.500 

1967.500 

1967.500 

30 

1964.250 

1964.250 

1964.250 

1962.083 

1961.991 

1962.668 

1962.783 

1965.850 

1965.650 

1965.850 

1965.850 

1965.850 

39 

1963. 500 

1963.500 

1963.500 

1961.049 

1960. H37 

1961.189 

1964.200 

1964.200 

1964.200 

1964,200 

1964.200 

1964,200 

40 

1963.000 

1963.000 

1963.000 

1982,000 

1961.034 

1982.000 

1963.000 

1963.000 

196 3.000 

1963.000 

1963.000 

1963.000 


ORIGINAi; PAGE 15 
OF POOR QUALITY 
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ITERATION - 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
TTEPATION = 
ITERATION = 
ITERATION ^ 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 

• ITERATION = 
ITERATION = 
ITERATION = 
ITERATION --- 
ITERATION ^ 
ITERATION = 
ITERATION = 
ITERATION ^ 
ITERATION - 
ITERATION = 
ITFRATIQN = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITERATION = 
ITEPATIOM = 
ITERAT ION = 
ITERATION - 


51 SUM= 

52 SUM= 

53 SUM= 

54 SUM= 

55 SUM= 

56 SUM= 

57 SUM^ 

50 SUiM= 

59 SUM= 

60 SUM= 

61 SUM= 

62 SUM= 

63 SUM= 

64 SUM= 

65 ■ 5UM= 

66 SUM= 

57 SUM= 

63 SUM= 

69 SUM= 

70 SUM= 

71 SUM- 

72 SUM- 

73 SUM- 

74 SUM- 

75 SUM- 

76 SUM- 

77 SUM- 

73 SUM- 

79 SUM- 

80 SUM- 

81 SUM- 

82 SUM- 

33 SUM- 

84 SUM- 

85 SUM- 

86 SUM- 

87 SUM- 

83 SUM- 

89 SUM- 

90 SUM- 

91 SUM- 

92 SUM- 

93 SUM- 

94 SUM- 

95 SUM- 

96 SUM- 

97 •' SUM- 

98 SUM- 

99 SUM- 

100 SUM- 


0.75354E 

02 

0.71660E 

02 

0.6B604E 

02 

0.65320E 

02 

0,62547E 

02 

0.59759E 

02 

0.57002E 

02 

0. 54649E 

02 

0.52126E 

02 

0.49842E 

02 

0.478C4E 

02 

0.45589E 

02 

0.436B4E 

02 

0.41854E 

02 

0.39955E 

02 

0.38292E 

02 

0.3671BE 

02 

Oa 35006E 

02 

0.33649E 

02 

0,32085E 

02 

0,308865 

02 

0,29443E 

02 

0,28252E 

02 

0„ 2T132E 

02 

0..25903E 

02 

0,248875 

02 

0,23816E 

02 

0,227345 

02 

Oo 21840E 

02 

0. 20906E 

02 

0.20005E 

02 

0,19153E 

02 

0, 184225 

02 

0,17635E 

02 

0. 16858E 

02 

0.16230E 

02 

0. 15522E 

02 

0. 14850E 

02 

0,14329E 

02 

0,135956 

02 

0 , 13 1 1 6 E 

02 

0. 12549E 

02 

0, 12022E 

02 

0, 11503E 

02 

0, 11082E 

02 

0, 1057BE 

02 

0. I0145E 

02 

0„97095E 

01 

0,93036E 

01 

0,89175E 

01 


45 


TFHPFRATUKE AT GRID POINTS 


J» 

I 

t 

2 

3 

• 

6 

5 

6 

7 

e 

9 

10 

11 

12 

1 

1961.000 

1961.000 

1961.000 

1968.999 

1968.609 

1969.437 

1961.000 

1961.000 

1961.000 

1961.000 

1961.000 

1961.000 

2 

1966.000 

1969.377 

1969.282 

1969.635 

1969. 778 

1970.541 

1971.797 

1972.674 

1998,500 

1998.500 

1998.500 

1998.500 

3 

1975.500 

1970.565 

1970.320 

1970.661 

1970.771 

1971.261 

1971.891 

1972.538 

1972.720 

1993.500 

1998.500 

1998.500 

A 

1970.979 

19 70.880 

1970.906 

1971.068 

1971.313 

1971.694 

1972.173 

1972.729 

1973.219 

1998.350 

1998.350 

1998,350 

5 

1975.200 

1971.928 

1971.289 

1971.660 

1971.681 

1972.007 

1972.609 

1972.869 

1973.360 

1973.519 

1998.200 

1998.200 

6 

1976.200 

1971.759 

1971.627 

1971.777 

1971.998 

1972.302 

1972.691 

1973.169 

1973.742 

1976.115 

199 8.150 

1998.150 

7 

1977.200 

19/2.153 

1972.052 

1972.121 

1972.288 

1972.545 

1972.899 

1973.357 

1971.939 

1976.556 

1998.100 

1998,100 

8 

1977.900 

1977.900 

1972.562 

1972.390 

1972.508 

1972.706 

1972.985 

1973.336 

1973.769 

1976.186 

1976.326 

1997.500 

9 

1978,600 

1978.600 

1972.821 

1972,699 

1972,717 

1972.844 

1973.063 

1973.370 

1973.757 

1976.219 

1976.488 

1996.900 

10 

1970,800 

1978.800 

1978.800 

1973.033 

1972. BOB 

1972.908 

1973.102 

1973.384 

1973.753 

1976.209 

1976.557 

1996.700 

11 

1979.000 

1979.000 

1979.000 

1973.002 

1972,817 

1972.911 

1973.086 

1973.346 

1973.689 

1976.117 

1976.522 

1996.500 

12 

1979.000 

1979.000 

1979.000 

1972.936 

1972,866 

1972.903 

1973.062 

1973.265 

1973.577 

1973.979 

1976.600 

1996.500 

13 

1979.000 

1979.000 

1979.000 

1979.000 

1973.107 

1972.878 

1972.961 

1973.134 

1973.396 

1973.765 

1976.220 

1996.500 

19 

1978.900 

1978.900 

1978.900 

1978.900 

1972.968 

1972.766 

1972.835 

1972.956 

1973.129 

1973.365 

1973.575 

1973,684 

15 

1978.800 

1978.800 

1978.800 

1978.800 

1972.953 

1972-829 

1972.617 

1972.880 

1973.009 

1973.191 

1973.409 

1973.576 

16 

1978.900 

1978.900 

1978.600 

1978.600 

1978.600 

1973.055 

1972.807 

1972.859 

1972.991 

1973.198 

1973.482 

1973.672 

17 

1978.000 

1978.000 

1978.000 

1978. 000 

1978.000 

1972.971 

1972.760 

1972.825 

1972.973 

1973.201 

1973.510 

1973.702 

18 

1976.900 

1976.900 

1976.900 

1976.900 

1976.900 

1972.824 

1972.662 

1972.766 

1972.9 3 2 

1973.166 

1973.478 

1973.669 

19 

1975.800 

1975.800 

1975. 800 

1975.800 

1975.800 

1972.699 

1972.583 

1972.746 

1972.909 

1973.079 

1973.299 

1973.470 

20 

1979.890 

1979.850 

1974.850 

1976.850 

1974.850 

1 972.896 

1972.806 

1973.015 

1973.125 

1973,077 

1972.710 

1972.817 

21 

1973,900 

1973.900 

1973.900 

1973.900 

1973.900 

1973.716 

1973.653 

1973.925 

1976.233 

1976.627 

1975.210 

1993.500 

22 

1973.050 

1973.050 

1973.050 

1973.050 

1973.050 

1973.050 

1976.827 

1974.873 

1975.176 

1975.587 

1976.073 

1992.500 

23 

1972.200 

1972.200 

1972.200 

1972.200 

1972.200 

1972.200 

1975.623 

1975.508 

1975.821 

1976.237 

1976.658 

1991.500 

29 

1971.700 

1971.700 

1971.700 

1971.700 

1971.700 

1971.700 

1975.765 

1975.867 

1976.182 

1976.586 

1976.932 

1990.050 

25 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1975.863 

1975.967 

1976.278 

1976.669 

1976.929 

1988.600 

26 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1975.728 

1975.852 

1976.161 

1975.495 

1976.635 

1986.900 

27 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1975.397 

1975.529 

1975.808 

1976.154 

197 6.220 

1985.200 

28 

1970.900 

1970.900 

1970.900 

1970.900 

1970.900 

1970.900 

1976.817 

1974.981 

1975.272 

1975.597 

1983.350 

1983.350 

29 

1970.600 

1970.600 

1970.600 

1970.600 

1970.600 

1973.882 

1973.888 

1974.164 

1976.670 

1974.684 

1981.500 

1981.500 

30 

1969.650 

1969.650 

1969.650 

1969.650 

1969.650 

1972.992 

1973.039 

1973,273 

1973.562 

1973.649 

1979.550 

1979.550 

31 

1968.700 

196 8.700 

1968.700 

1968.700 

1968.700 

1972.068 

1972.165 

1972.330 

1972.670 

1977.600 

1977.600 

19T7.600 

32 

1968.350 

1968.350 

1968.350 

1968,350 

1968.350 

1971.282 

1971.391 

1971.592 

1971.758 

1975.800 

19TS.800 

1975.800 

33 

1968.000 

1968.000 

1968.000 

1968.000 

1968.000 

1970.461 

1970. 593 

1970.776 

19T0.863 

1974.000 

1976.000 

1974.000 

39 

1967.250 

1957.250 

1967,250 

1967.250 

1969.546 

1969.568 

1969.733 

1969.905 

1972.500 

1972.500 

1972.500 

1972.500 

35 

1966.500 

1V56.500 

1966.500 

1966,500 

1968.696 

1968.752 

1968.920 

1969.009 

1970.500 

1970.500 

1970.500 

1970.500 

36 

1965.750 

1965,750 

1965.750 

1965.750 

1957.90 3 

1966.015 

1968.267 

1968.281 

1969.000 

1969.000 

1969.000 

1969.000 

37 

1965.000 

1965.000 

1965.000 

1965.000 

1966.843 

1966.975 

1967.066 

1967.500 

1967.500 

1967.500 

1967.500 

1967.500 

38 

1969.250 

1966.250 

1966.250 

1965.886 

1965.890 

1965.966 

1965.984 

1965.850 

1965.850 

1965.850 

1965.850 

1965.850 

39 

1963. 50 0 

1963. 500 

1963. 500 

1965.565 

196 5.60 5 

1965.618 

1964.200 

1966.200 

1966.200 

1964.200 

1964.200 

1964.200 

90 

1963.000 

1963.000 

1963.000 

1982.000 

1965.501 

1982,000 

1963.000 

1963.000 

1963.000 

1963.000 

1963.000 

1963.000 
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ITERATIGM 
ITERATION 
I TERATION 
ITERATION 
ITEPATIOM 
ITERATION 
ITERATION 
ITERATION 
ITERAT ION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERAT TOM 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
I TERATION 
ITERAT ION 
ITERATION 
ITERATION 
ITERAT ION 
I TERATION 
ITERAT ION 


= L0I 

SUM= 0,a587BE 

01 

= 102 

SUM= Oo8l821E 

01 

= 103 

SUN= 0e7Ba01E 

01 

= 10't 

SUM= 0o75410E 

01 

=105 

SUM= 0.72450E 

01 

= 106 

SUM= 0.690 65 E 

01 

= 107 

SUM= 0o66719E 

01 

= 108 

SUH= 0„63550E 

01 

= 109 

SUH= 0«61287E 

01 

= 110 

SUM= 0»58460E 

01 

= 111 

SU,M= 0.56228E 

01 

= l 12 

5UM= 0.53530E 

01 

= 113 

SUM= 0.51897E 

01 

= 114 

SUM= 0.49356E 

01 

= 115 

SUM= 0.47412E 

01 

= 116 

SUH= 0.45443E 

01 

= 117 

5UH= 0.A3495E 

01 

= 118 

SUM= 0o4L741E 

01 

= 119 

SUM= 0,40085E 

01 

= 120 

SUM= 0.38259E 

01 

= 121 

SUM= 0,36889E 

01 

= 12? 

SUM= Oo35259E 

01 

= 123 

SUM= 0„33984E 

01 

= 124 

SU,M= 0a32260E 

01 

= 125 

SIJM= 0.30942E 

01 

= 126 

SUM= 0.29916E 

01 

= 127 

SU.M= 0,28456E 

01 

= 128 

SUM= 0.27564E 

01 

= 129 

SUM= 0o26040E 

01 

= 130 

SUM= 0«25450E 

01 

= 131 

SUM= 0.23825E 

01 

=132 

SUM= 0<,23447E 

01 

= 133 

SUH= 0.220026 

01 

= 134 

5UM= 0.21550E 

01 

= 135 

SUH= 0.20353E 

01 

= 136 

SUN= 0.19933E 

01 

= 137 

SUM= 0.18697E 

01 

= 13 8 

SUM= 0.18322E 

01 

= 139 

SUiM= 0.17353'E 

01 

= 140 

SUM= 0.L6836E 

01 

= 141 

SIJM= 0. 15795E 

01 

= 142 

SUM= 0,153586 

01 

= 143 

• SUM= 0. 14590 E 

01 

= 144 

SUM= 0.13990E 

01 

= 145 

SUM= 0.13506E 

01 

= 146 

SUN= 0,127736 

01 

= 147 

SUH= 0,12664E 

01 

= 148 

SUiM= 0-11726E 

01 

= 149 

SUM= 0.11372E 

01 

= 150 

SUN= 0.10099E 

01 
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TFMPEKATURE AT GRID POINTS 


1 

1 

2 

3 

9 

5 

6 

7 

8 

9 

10 

u 

12 

1 

1961.000 

1961.000 

1961.000. 

1979.976 

1979. 159 

1979. 833 

■961.000 

1961.00 0 

1961. 000 

1951.000 

1961.000 

1961.000 

2 

1966.000 

1975.903 

1975.392 

1975.583 

1975.901 

1976.566 

1977.691 

1978.361 

1998.500 

1993.500 

1998. 306 

1998.500 

3 

1975.500 

1976.879 

1976.789 

1976.929 

1977.218 

1977.656 

1978.212 

1978.791 

1978.908 

1998.500 

1998.500 

1998.500 

4 

1977.569 

1977.591 

1977.612 

1977.777 

1978.035 

1978.378 

1978.790 

1979.298 

1979.610 

1998.350 

1998.350 

1998,350 

5 

1,975.200 

1978.079 

1978. 160 

197P.363 

1978.617 

1978.927 

1979.286 

1979.683 

1980.110 

1960.198 

1998.200 

1998.200 

6 

1976.200 

1978.611 

1978.697 

1978.899 

1979.093 

1979.309 

1979,725 

1980. 115 

1980.560 

1960.306 

1998. 150 

1998.150 

7 

1977.200 

1979.132 

1979. 123 

1979.202 

1979.991 

1979.751 

1980.065 

■ 1930.936 

1980.879 

1981.297 

1993.100 

1993.100 

U 

1977.900 

1 977.900 

1979.587 

1979.623 

1979.799 

1980.019 

1980.289 

1980.599 

1980.995 

1981.323 

1981.390 

1997.500 

9 

1978.600 

1978.600 

1979.930 

1979.919 

1930.039 

1980.213 

1980.937 

1980.706 

1981.011 

1931.356 

1981.512 

1996.900 

10 

1978.800 

1978.800 

1978.600 

1980.163 

1980.177 

1980.327 

1960.522 

1980.762 

1981.091 

1981.360 

1981.569 

1996.700 

11 

1979.000 

1979.000 

1979.000 

1980.298 

1980.238 

1980.375 

1980. 598 

1980.762 

1931.015 

1981.306 

1981. 559 

1996.500 

12 

1979.000 

1 979.000 

1979.000 

1980.326 

1980.295 

1980.386 

1980.527 

1980.719 

1930.999 

1991.215 

1931.973 

1996.500 

13 

1979.000 

1979.000 

1979.000 

1979.000 

19R0.312 

1980.305 

1980.935 

1980.605 

1930.819 

1981.065 

1931.353 

1996.500 

lA 

1978.900 

1 978.900 

1978.900 

1978.900 

1980.189 

1980. 166 

1980.293 

1980.997 

1980.628 

1980.332 

1931.095 

1931.069 

15 

1978. SOO 

1978.800 

1976.800 

1978.800 

1980.088 

1980.069 

1930.158 

1980.299 

1980.963 

1980.660 

1980.879 

1980.959 

16 

1978.900 

I 978.900 

1978.900 

1978.900 

1978.900 

1979,976 

*979.981 

1960.116 

19B0.295 

1980.510 

1980.765 

1980.862 

17 

1978.000 

1978.000 

1978.000 , 

1978.000 

1978.000 

1979.730 

1979.733 

1979.682 

1980.075 

1980.306 

1980.570 

1980. 67S 

18 

1976.900 

1976.900 

1976.900 

1976.900 

1976.900 

1979.395 

1979.903 

1979.582 

1979,799 

1980. P91 

1930.325 

1980,927 

19 

1975.800 

1975.800 

1975.800 

1975.800 

1975.800 

1979.022 

1979.029 

1979.252 

1979.9 79 

1979.718 

1979.979 

1980.067 

20 

1979.850 

1979.850 

1979.850 

1979.050 

1979.850 

1978.732 

1978.738 

1978.991 

1979.213 

1979.386 

1979.958 

1979, ,501 

21 

1973.900 

1973.900 

1973.900 

1973. 9C0 

1973.900 

1978.658 

1978.653 

1978 .937 

*.97«!.?39 

’.979,56?. 

l'>79.959 

1993.500 

22 

1973.050 

1973.050 

1973.050 

1973.050 

1973.050 

1973.050 

1978.673 

1978.885 

1979.195 

1979.598 

1979.915 

1992.500 

23 

1972.200 

1972.200 

1972.200 

1972.200 

1972.200 

1972.200 

1978.500 

1978.709 

1979.035 

1979.903 

1979.739 

1991.500 

2A 

1971.700 

1971.700 

1971.700 

1971.700 

1971.700 

1971.700 

1978,219 

1978.923 

1978.750 

1979.116 

1979.399 

1 990.050 

25 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1977.810 

1978.020 

1978.399 

1978.705 

1978. 91y 

1936.600 

26 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1977.311 

1977.503 

1977.805 

1978.199 

1978,256 

1986.900 

27 

1971.200 

1971,200 

1971.200 

1971.200 

1971.200 

1971.200 

1976.651 

1976.899 

1977.130 

1977.960 

1977.508 

1985.200 

28 

1970,900 

1 970.900 

1970.900 

1970.900 

1970.900 

1970.900 

1975.802 

1976.022 

1976.309 

1976.602 

1983.350 

1983.350 

29 

1970.600 

1970.600 

1970.600 

L970.600 

1970.600 

1979.790 

1979.762 

1975.029 

1975.312 

1975.509 

1981.500 

1981.500 

30 

1969.650 

1 969. 650 

1969.650 

1969.650 

1969.650 

1973,633 

1973,700 

1973.937 

1979.190 

1979.235 

1979.550 

1979.550 

31 

1968.700 

1968.700 

1968.700 

1968.700 

196 8.700 

1972.529 

1972.626 

1972.300 

1972.90- 

1977.600 

1977.600 

1977.600 

32 

1968.350 

1968.350 

1968.350 

196 8.350 

1960.350 

1971.693 

1971,772 . 

1971. 972 

1972.126 

1975.800 

1975.800 

1975. SOO 

33 

1 568.000 

1968.000 

1968.000 

1968.000 

1968.000 

1970. 758 

1970.908 

1971.088 

1971.152 

1979.000 

1979.000 

1979.000 

3A 

1967.250 

1967.250 

1967.250 

1967.250 

1969.322 

1969.833 

1970.005 

1970.166 

1972.500 

1972.500 

1972.500 

1972.500 

35 

1966.500 

1966.500 

1966.500 

1966. 500 

1968.905 

1968.971 

1969.139 

1969.216 

1970.500 

1970.500 

1970.500 

1970.500 

36 

1965. 750 

1965. 750 

1965.750 

1965.750 

1968.063 

1968. 181 

1968.398 

1968.911 

1969.000 

1969.000 

1969.000 

1969.000 

37 

1965.000 

1965.000 

1965.000 

1965.000 

1967.063 

1967.212 

1967.288 

1967.500 

1967.500 

1967.500 

196 7.500 

1967. 500 

38 

1969.250 

1969.250 

1969.250 

1966. 10 6 

1966. U6 

1966.131 

1966, 167 

1965.850 

1965.850 

1965.850 

1965.850 

1965.850 

39 

1963.500 

1963.500 

1963.500 

1965. S21 

1965.079 

1965.870 

1969.200 

1969.200 

1969.200 

1969.200 

1969.200 

1964.200 

60 

1963.000 

1963.000 

1963.000 

1982.000 

1966.759 

1982.000 

1963.000 

1963.000 

1963.000 

1963.000 

1963.000 

1963.000 
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ITERATION 

ITERATION 

ITERATION 

ITERATION 

ITERATION 

ITERATION 

ITERATtON 

ITERATION 

ITERATION 

ITERATION 

iteration 

ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITER AT ION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
HERAT ION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERATION 
ITERAT K3N 
ITERATTON 
ITERATION 
ITERATION 
ITERATION 
ITERATION 


=151 

SUM= 0.10474E 

01 

= 152 

SUM= 0.10058E 

01 

= 153 

SUM= 0.97168E 

00 

= 154 

SUM= 0.91090E 

00 

= 1 55 

SUM= 0.90131E 

00 

= 156 

SUM= 0.84379E 

00 

= 157 

SUM= 0,82349E 

00 

= l5fT 

SUM= 0.77830E 

00 

= 159 

SUM= 0.74847E 

00 

= 160 

SUM= 0.71357E 

00 

= 161 

SUM= 0.60149E 

00 

= 162 

SUM= 0.6657LE 

00 

= 163 

SUM= 0.62334E 

00 

= 164 

SUM= 0.61099E 

00 

= 165 

SIJM= 0.58113E 

00 

= 166 

SUM= 0.55107E 

00 

= 167 

SUM= 0.54886E 

00 

= 168 

SUM= 0,51087E 

00 

= 169 

SUM= 0o50309E 

00 

= 170 

SUM= 0.47725E 

00 

= 171 

SUM= 0,46802E 

00 

=172 

SUiM= 0.43639E 

00 

=173 

SUM= 0-.43087E 

00 

= 174 

SUM= 0.40216E 

00 

= 175 

SUM= 0.38674E 

00 

= 176 

SUM= 0.37460E 

00 

= 177 

SUM= 0.35311E 

00 

=170 

SUM= 0.342 74 E 

00 

=179 

SON= 0.33164E 

00 

= 100 

SUM= 0.3T819E 

00 

= 181 

SUN= 0.30598E 

00 

= 182 

SUM= 0.28902E 

00 

= 183 

SUH= 0.28203E 

00 

= 1 84 

SUM= 0.26825E 

00 

= 105 

SUM= 0.25814E 

00 

= 136 

SUM= 0.25256E 

00 

= 1 87 

SUM= 0.24203E 

00 

= 103 

SUM= 0.22884E 

00 

=189 

SUM= 0.21402E 

00 

= 190 

SUN= 0.21737E 

00 

= I9i 

SUM= 0.20035E 

00 

= 192 

5UN= 0.20713E 

00 

= 193 

SUM= 0.19532E 

00 

= 194 

SUM= 0.18689E 

00 

= 1 95 

SUM= 0.18265E 

00 

= 196 

SUM= 0.16749E 

00 

= 197 

SUM= 0-16794E 

00 

= 198 

S'JM= 0.15877E 

00 

= 199 

SUM= 0.1492BE 

00 

=^ 00 

5UM= 0. 14393E 

00 
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TEP*PEKATURE AT GRID POINTS 


J» 

1 

1 

2 

3 . 

4 

5 

6 

7 

8 

9 

10 

U 

12 

1 

1961.000 

1961.000 

1961.000 

1976.419 

1976. U3 

1976.750 

1961.000 

1961.00 0 

1961.000 

1961.000 

1961.000 

1961.000 

2 

1966.000 

1977.541 

1977.560 

1977.765 

1978,075 

1970.706 

1979.716 

1980.380 

1998.500 

1998.500 

1998.500 

1998.500 

3 

1975.500 

1979.113 

1979.082 

1979.224 

1979.505 

1979.926 

1980.456 

1981.010 

1981.103 

1998.500 

1998. 500 

1998.500 

A 

1579.900 

1 979.905 

1979.991 

1900.165 

1980,420 

1980.750 

1981.133 

1981.560 

1981.877 

1998.350 

1998.350 

1998.350 

5 

1975.200 

1980.433 

1980.595 

1980. 816 

191)1. 076 

1981.379 

1981.721 

1982.097 

1982.499 

1982.559 

1998.200 

1998.200 

6 

1976.200 

1981.034 

1981.130 

1981.351 

1981.604 

1981.889 

1902.211 

1982.571 

1982.969 

1983.169 

1998. 150 

1998. ISO 

7 

1977.200 

1981.597 

1981.619 

1981.809 

1982.034 

1982.295 

1982.594 

1982.934 

1983.319 

1983.673 

1998.100 

1998.100 

8 

1977.900 

1977.900 

1982.061 

1962.169 

1982.365 

1982.593 

1982.656 

1983.151 

1983.479 

1983.836 

1983.875 

1997.500 

9 

1978.600 

1978.600 

1982.433 

1982.459 

1982.613 

1982.804 

1983.027 

1933.280 

1933.561 

1983.865 

1983.962 

1996.900 

10 

1978.800 

1978.800 

1978.800 

1982.659 

1982.757 

1982.926 

1983.122 

1983.345 

1983.592 

1983.863 

1984.023 

1996,700 

11 

1979.000 

1979.000 

1979.000 

1982.779 

1982.830 

1962.981 

1983.152 

1983.348 

1983.570 

1983.812 

1984.003 

1996.500 

12 

1979.000 

1979.000 

1979.000 

1982.899 

1982.888 

1982.989 

1963.129 

1963.302 

1983.502 

1983.726 

1983.926 

1996.500 

13 

1579.000 

1979.000 

1979.000 

1979.000 

1982.807 

1982.876 

1983.020 

1983.186 

1983.375 

1983,587 

1983.814 

1996.500 

H 

1978.900 

1978.900 

1978.900 

1978.900 

1982. 674 

1982.712 

1982.856 

1983.019 

1983.201 

1983.400 

1983.611 

1983.631 

15 

1978.800 

1978. BOO 

1978. 800 

1978.800 

198 2.524 

1982.536 

1982.667 

1962.823 

1983.004 

1983.205 

1983.423 

1983.475 

16 

1970.400 

1978.400 

1978.400 

1978.400 

1978.400 

1982. 320 

1982.412 

198 2.57 5 

1982.767 

1982.983 

1983.238 

1983.295 

17 

1978.000 

1978.000 

1978.000 

1978.000 

.1978.000 

1982.004 

1982.079 

1982.254 

1992.460 

1982.692 

1982.952 

1983.021 

18 

1976.900 

1976.900 

1976.900 

1976.900 

1976.900 

1981.589 

1981.653 

1981.857 

1982.083 

1982.333 

1982.608 

1982.679 

15 

1975.800 

1975. 800 

1975.800 

1975. 800 

1975.800 

1981.116 

1901. 163 

1981.404 

1981.652 

1981.911 

1982.188 

1982.249 

20 

1974.850 

1974. 850 

1974.850 

1974.850 

1974..85U 

1980.649 

1980.667 

1980.955 

1931.212 

:9Q 1.456 

1981-666 

1981. 689 

21 

1973.900 

1973. 90 0 

1973.900 

1973.900 

1973.900 

1980.273 

1980.287 

1980.574 

1980.668 

1981,176 

1981.505 

1993.500 

22 

1973.050 

1973.050 

1973.050 

1973.050 

1973.050 

1973.050 

1979.923 

1980.187 

1980.499 

1980.832 

1981.160 

1992.500 

23 

1972.200 

1972.200 

1972.200 

1972.200 

1972.200 

1972.200 

1979.489 

1979.739 

1980.068 

1980.419 

1980.721 

1991.500 

24 

1971.700 

1971.700 

1971.700 

1971.700 

1971. 700 

1971,700 

1973.988 

1979.230 

1979.563 

1979.918 

1980.175 

1990.050 

25 

1971.200 

1971.200 

1971.200 

1971,200 

1971.200 

1971.200 

1978.422 

1978.658 

1978.985 

1979.340 

1979.540 

1988.600 

26 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1977.792 

1978.004 

1978.307 

1978.643 

1978.746 

1986.900 

27 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1977.017 

1977.228 

1977.518 

1977.847 

1977.892 

1985.200 

28 

1970.900 

1970.900 

1970.900 ' 

1970.900 

1970.900 

1970.900 

1976.075 

1976.310 

1976.598 

1976.804 

1983.350 

1983. 35C 

25 

1970.600 

1970.600 

1970.600 

1970.600 

1970.600 

1974.963 

1974.986 

1975.248 

1975.532 

1975.720 

1981,500 

1981.500 


,969,650 

1969.650 

1969.650 

1969.650 

1969.650 

1973.791 

1973.861 

1974.101 

1974.352 

1974.445 

1979.550 

1979.550 

■Ji 

1968.700 

196B.TOO 

1968.700 

1968.700 

1968.700 

1972.636 

1972.73 T 

1972.909 

1973.001 

1977.600 

1977.600 

1977.600 

32 

1968.350 

1968.350 

1968.350 

1968.350 

1968.350 

1971.723 

1971.855 

1972.054 

1972.207 

1975.800 

1975.800 

19f5.80C 

33 

1968.000 

1968.000 

1968.000 

1968.000 

1968.000 

1970.818 

1970.972 

1971.153 

1971.215 

1974.000 

1974.000 

1974.000 

34 

1967.250 

1967.250 

1967.250 

1967.250 

1969.879 

1969.690 

1970.060 

1970.219 

1972. 500 

1972.500 

1972.500 

1972.500 

35 

1966.500 

1966. 500 

1966.500 

1966.500 

1968.943 

1969.011 

■ 1969.173 

1969.256 

1970.500 

1970.500 

1970.500 

1970.500 

36 

1965.750 

1965.750 

1965.750 

1965.750 

1968.087 

1968.206 

1968.424 

1969.436 

1969.000 

1969.000 

1969.000 

1969.000 

37 

1965.000 

1965.000 

1965.000 

1965.000 

196 7.080 

1967.230 

1967.307 

1967.500 

1967.500 

1967.500 

1967.500 

1957.500 

36 

1964.250 

1964.250 

1964.250 

1966.108 

1966.119 

1986.136 

1966.174 

1965.850 

1965.850 

1965.850 

1965.850 

1965.850 

39 

1963.500 

1963. 500 

1963.500 

1965.825 

1965.880 

1965.975 

1964.200 

1964.200 

1964.200 

1964.200 

1964.200 

1964.200 

40 

1963.000 

1963.000 

1963.000 

1982.000 

1965.759 

1982.000 

1963.000 

1963.000 

1963.000 

1963.000 

1963.000 

1963.000 


50 


nP POOR ODAlirS 


, HERAT ION 

=201 

ITERATION 

=202 

I TERATIOM 

= 203 

ITERATION 

=2 04 

ITERATION 

=205 

ITERATION 

=2 06 

ITERATION 

=2 0 7 

I TER AT ION 

=2 0 0 

ITERATION 

= 2 09 

ITERAT ION 

= 210 

ITERATION 

= 211 

ITERATION 

= 212 

ITERATION 

=213 

ITERATION 

= 214 

ITERATION 

= 215 

ITERATION 

= 216 

ITERATION 

=217 

ITERATION 

= 218 

I TERATION 

= 219 

ITERAT ION 

=220 

ITERATION 

=2 21 

ITERATION 

= 222 

ITERATION 

=223 

ITERATION 

=2 24 

ITERAT ION 

=225 

ITERATION 

= 2 26 

ITERATION 

= 22 7 

ITERATION 

= 228 

ITERATION 

=2 29 

ITERATION 

=2 30 

ITERATION 

= 23L 

ITERATION 

=2 32 

ITERATION 

=2 33 

ITERATION 

= 234 

ITERATION 

= 235 

ITERATION 

=2 3 6 

ITERATION 

=23 7 

ITERATION 

=238 

ITERATION 

=239 

ITERATION 

=240 

ITERATION 

= 241 

ir ERAT lOM 

= 292 

ITERATION 

= 2 93 

ITERATION 

=2 49 

HERAT ION 

=2 45 

ITERATION 

= 2 46 

ITERATION 

=247 

ITERATION 

=248 

ITERATION 

=249 

ITERATION 

=2 50 


SUM= 0.13653E 00 
SUM= 0,13276E 00 
SUH= 0.12584E 00 
SL;M= 0. UHOE 00 
SUM= Ow 12049E 00 
SUM= 0.12342E 00 
SUM--^ 0*11000': 00 
SUM- 0*10766S 00 
SUM= 0. 10171E 00 
SUM= 0.1016LE no 
SUM= 0.91633E-01 
SUM= 0.90S95E-01 
SUM= 0.848565-01 
SUM= 0.87407E-01 
SUM= 0. 810525-01 
SUM= 0.84283E~01 
SUM= 0,80870E-01 
SUM= 0.73778E-01 
SUM= 0*741465-01 
SUM= 0.708555-01 
SUM= 0. 663305-01 
SUM= 0.64120E-01 
S1JM= 0.619725-01 
SUM= 0.61971E-01 
SUiM= 0.656745-01 
SUM= 0.57977E-01 
SUN5= 0.56381E-01 
SUM= 0.52466E-01 
SUM= 0.519625-01 
SUM= 0.535385-01 
.SUi^= 0.54S15E-01 
SOM= 0.56712E-01 
SUM= 0.5641CE-01 
SlJM= 0.53055E-01 
SUM- 0.48655E-01 
SUM= 0.43210E-01 
5UM= 0.46932E-01 
SUM= 0.535a4£-01 
SUM= 0-46352E-01 
SUM= 0.43471E-01 
SUM= 0.914235-01 
SUM= i) .4 38 23 5-01 
SUM= 0.37151=-01 
SUM= 0.36S63E-01 
SUM= 0.906205-01 
5UM= 0.929715-01 
SUM= 0.990845-01 
SUM- 0.37393E-01 
SUM= 0.375465-01 
SUM= 0.373005-01 
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TF*1PFRATURE AT GRin POINTS 


J* 

1 

1 

2 

3 

• 

4 

5 

6 

7 

e 

9 

10 

11 

12 

1 

1961.000 

1961.000 

1961 .000 

1977.104 

1976.01 r 

1977.423 

J961.000 

1961.000 

1961.000 

1961.000 

1961.000 

1961.000 

2 

1966.000 

1978.295 

1978.325 

1978.535 

1978.841 

1979.458 

1900.444 

1981.089 

1998.500 

199R.500 

1998.500 

1998,500 

3 

1975.500 

1979.904 

1979.887 

1980.031 

1980.309 

1980.722 

1981.243 

1981.786 

1981. 870 

1998.500 

1998.500 

1998.500 

4 

1980.716 

1980. 730 

1980,823 

1980,999 

190 U252 

1981.577 

1981.958 

1982.367 

1982. 66B 

1990.350 

1998.350 

1998.350 

5 

1975.200 

1981.255 

1981.446 

1981.673 

1931.934 

1982.235 

1982.573 

1982.939 

1983.335 

1903.308 

1998.200 

1998.200 

b 

1976.200 

1981.802 

1981.999 

1982.226 

1982.478 

1982.762 

1963.079 

1983.426 

1983. 806 

1983.990 

1998.150 

1998.150 

7 

1977.200 

1982,456 

1982.491 

1982.691 

1982,921 

1983.181 

1983.474 

1963.802 

1984. 169 

1904,499 

1990. 100 

1998.100 

8 

1977.900 

197 7.900 

1982.922 

1983.055 

1983.258 

1983.490 

1983.749 

1984.039 

1984.357 

1904.710 

I9B4.743 

1997.500 

9 

1978.600 

1978.600 

1983,305 

1983.344 

1983.508 

1983.703 

1983.925 

1984.175 

1964.447 

1984.739 

1904.841 

1996.900 

10 

1978. 800 

1978.800 

1978.800 

1983.527 

1983.652 

1983.025 

1984.020 

1984.239 

1984.475 

1964.730 

1984.873 

1996.700 

11 

1979.000 

1 979,000 

1979.000 

1983.651 

1983.723 

1903. 8/7 

1904.049 

1984.240 

1984.449 

1984.674 

1984.847 

1996.500 

12 

1979.000 

1979.000 

1979.000 

1983.778 

1983.775 

1983.880 

1904.020 

1984.188 

1984.377 

1904.585 

1984.765 

1996.500 

13 

1979.000 

1 979.000 

1979,000 

1979.000 

1983.655 

1983.751 

1983.900 

1984.068 

1984.251 

1984.448 

1984.654 

1996.500 

14 

1978.900 

1978.900 

1978.900 

1978.900 

1983.523 

1983.579 

1983.732 

1983.897 

1984.080 

1984.277 

1984. 4H4 

1984.497 

15 

1978.800 

1978.800 

1978.800 

1978.600 

1983.355 

1983.377 

1983.519 

1983.684 

1903.869 

1984.073 

1904.293 

1984.335 

16 

1978.400 

1978.400 

1978.400 

1978.400 

1976.400 

1983.117 

1983.236 

1983.407 

1983.603 

1983.823 

1984. 064 

1984.120 

17 

1978.000 

1978.000 

1978.000 

1978.000 

1978.000 

1982.771 

1982.871 

1983.054 

1983,264 

1983.498 

1983. 753‘ 

1983. B12 

18 

1976.900 

1976.900 

1976,900 

1976.900 

1976.900 

1982.326 

1982.409 

1982.620 

19B2.852 

1983.103 

1983.375 

1963.436 

19 

1975.800 

1975.800 

1975.800 

1975.800 

1975. 800 

1981.817 

1981.877 

198 2.124 

19B2.30O 

1962.645 

1982.926 

1982.979 

20 

1974.350 

1974.850 

1974-8.50 

5974.850 

l97i,H50 

1981.791 

1081.333 

1981.609 

1981. PRO 

1982.147 

196?, 403 

5982.419 

21 

1973.900 

1973. 900 

1973.900 

1973. 900 

1973.900 

1980.810 

1980. 830 

1901.118 

1961,411 

1981.711 

1982.020 

1993.500 

22 

1973.050 

1973.060 

1973.050 

1973.050 

1973.050 

1973.050 

1980.336 

1980.618 

1980.932 

1981,259 

1981.575 

1992.500 

23 

1972.200 

1972.200 

1972.200 

1972.200 

1972.200 

1972.200 

1979.813 

1980.075 

1980.407 

1980.755 

1981.049 

1991.500 

24 

1971.700 

1971.700 

1971.700 

1971.700 

1971.700 

1971.700 

1979.240 

1979.493 

1979.827 

1980. 178 

1980.428 

1990.050 

25 

1971.200 

1971.200 

1971.200 

1971. ZOO 

1971,200 

1971.200 

1978.618 

1978.863 

1979.190 

1979.542 

1979.737 

1988.600 

26 

1971.200 

1971.200 

1971.200 

1971.200 

1971,200 

1971.200 

1977.943 

1978.161 

1978,466 

1978.800 

1978.902 

1986.900 

27 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1977.133 

1977.347 

1977,639 

1977.967 

1978 , 0 10 

1985.200 

28 

1970.900 

1970.900 

1970.900 

1970.900 

1970.900 

1970.900 

1976.156 

1976.395 

1976.685 

1976.969 

lv83.350 

1983.350 

29 

1970.600 

1970.600 

1970.600 

1970.600 

1970.600 

1975.008 

1975.035 

1975.303 

1975,591 

1975.777 

1981.500 

1961.500 

30 

1969.650 

1969.650 

1969.650 

1960.650 

1969.650 

1973.827 

1973.898 

1974.140 

1974.392 

1974.485 

1979,550 

1979.550 

31 

1968.700 

1968.700 

1968.700 

1968.700 

1968.700 

1972.657 

1972.758 

, 1972.932 

1975.022 

197 7.600 

1977.600 

1977.600 

32 

1968.350 

1968.350 

1968.350 

1968.350 

1960.350 

1971.738 

1971.870 

1972.071 

1972.224 

1975,800 

1975.800 

1975.800 

33 

1968.000 

1968.000 

1968. QUO 

1968.000 

1968.000 

1970.826 

1970.982 

1971.162 

1971.224 

1974.000 

1974,000 

1974.000 

34 

1967.250 

1967.250 

1967.250 

1967.250 

1969.881 

1969.692 

1970.062 

1970.222 

1972.500 

1972,500 

1972.500 

1972.500 

35 

1966.500 

1966.500 

1966.500 

1966. 500 

1960.944 

1969.010 

1969.174 

1969.257 

1970.500 

1970.500 

1970.500 

19TO.SOO 

36 

1965. 750 

1965.750 

1965.750 

1965.750 

1960.093 

1966.213 

1968.428 

1968.438 

1969.000 

1969.000 

1969.000 

1969.000 

37 

1965.000 

1965,000 

1965.000 

1965.000 

1967.005 

1967.237 

1967.313 

1967.500 

1967.500 

1967.500 

1967.500 

1967.500 

38 

1964.250 

1964.250 

1964.250 

1966. 113 

1966.124 

1966.142 

1966.177 

1965.850 

1965.850 

1965.850 

1965.850 

1965.850 

39 

1963.500 

1963.500 

1963.500 

1965.833 

1965.887 

1965.882 

1964.200 

1964.200 

1964.200 

1964.200 

1964.200 

1964.200 

40 

1963.000 

1963.000 

1963.000 

1982.000 

1965.766 

1982.000 

1963.000 

1963.000 

1963.000 

1963.000 

1963.000 

1963.000 


original PAG^ 
OR POOR QUALH^ 


ITERATION =251 
ITERATION =252 
ITERATION =253 
ITERATION =254 
iteration =255 
ITERATION =256 
ITERATION =257 
ITERATION =253 
ITERATION =259 
ITERATION =260 
ITERATION =261 
ITERATION =26? 
ITERATION =263 
ITERATION =264 
ITERATION =265 
ITERATION =266 
ITERATION =267 
ITERATION =263 
ITERATION =269 
ITERATION =270 
ITERATION =271 
ITERATION =272 
ITERATION =273 
ITERATION =274 
ITERATION =275 
ITERATION =276 
ITERATION =277 
ITERATION =278 


SUM= 

0, 

.38970E- 

■01 

SUM= 

0. 

.35360E- 

•01 

SUM= 

0. 

. 34429E- 

■01 

SUM= 

0. 

.411A0E- 

■01 

SUM= 

0. 

. 332340- 

■01 

SUM= 

0, 

,33365E- 

01 

SUM= 

0, 

.32531E- 

•01 

SUM= 

0. 

►26938E- 

•01 

SU.M= 

0, 

.265485- 

01 

SUM= 

0. 

.30519E- 

01 

SUM = 

0, 

, 25117E- 

01 

SUM = 

0, 

•27695E- 

01 

SUM= 

0, 

.28535E- 

01 

SUM= 

0, 

.251L5E- 

01 

SUM= 

0 , 

.28753E- 

01 

SIJM,= 

0 . 

.26751E- 

01 

5UM= 

0«233d6E- 

01 

SUM= 

0, 

. 25677 E- 

01 

SOM = 

0. 

.23795E- 

01 

SUiM= 

0, 

. 2498.5 E- 

01 

SUM= 

0, 

,21331E- 

01 

SUM= 

0<, 

■22153E- 

01 

SUM = 

0- 

.21594 E- 

01 

SUM= 

0, 

,23230E- 

01 

SUM = 

0, 

.25253E- 

01 

SUN= 

0. 

.27588E- 

01 

SUM= 

0, 

24965E- 

01 

SUM= 

Oo 

.19849E- 

01 
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THE FINAL ITtKATIVE 


TEMPERATURE AT GRID POINTS 


J= 

I 

1 

2 

3 

A 

5 

6 

7 

8 

9 

10 

11 

12 

1 ' 

1961.000 

1961.000 

1961,000 

1977. 26A 

1976.919 

1977. 580 

1961.000 

1961.000 

1961.000 

1961.000 

1961.000 

1961.000 

2 

196 6. 000 

1978. A69 

1978.503 

1978.715 

1979.020 

1979. 63A 

1980.615 

19B1.255 

1998.500 

1998.500 

1996. 500 

1998.500 

3 

1975.500 

1980.090 

1980.077 

1980,219 

1980, A98 

1960.910 

1981. A27 

1981.971 

1902.055 

1998.500 

1990.500 

1998.500 

A 

1980. 912 

1980.928 

1981.021 

1981.197 

1901. A51 

1981. 77A 

1962.151 

1982.558 

1982.856 

1998.350 

1998.350 

1998.350 

5 

1975.200 

1981. A51 

1961. 6A9 

1981.876 

1982.138 

1982. A37 

1982.773 

1983.140 

1983.532 

1983.583 

1998.2M 

1998.200 

t 

1976.200 

1982.083 

19d2.20A 

1982. A32 

1982.687 

1982.969 

19B3.28A 

1983.629 

1984.007 

1984. 188 

1998, 150 

1 998.150 

7 

1977. 200 

1982. 662 

1982.698 

1982.900 

1983.132 

1983.392 

1983.603 

1984.009 

1984.371 

1.984,695 

199$v 100 

1998.100 

S 

1977.900 

1977.900 

1983.129 

1963.267 

1983. ATI 

1983.703 

1903.962 

1984.250 

1984.565 

) 984,914 

f 92 4. 944 

1997.500 

s 

1978.600 

1978.600 

1983.510 

1983.552 

1983.720 

1983.916 

1984.139 

1984.385 

1984, 6 

19.J4.940 

198 5.041 

1 9«6.900 

10 

1978. 800 

1978.800 

1978. 800 

1983.730 

1983. 063 

198A.037 

19BA.232 

1984.447 

1984, 6Si 

1964.932 

1965.071 

1996.700 

11 

1979, 000 

1979.000 

1979.000 

1983.857 

1983.933 

198A.0S9 

1984. 259 

1984.449 

1984.655 

1984.876 

1985. 043 

1996.500 

12 

1979.000 

1979.000 

1979.000 

1983.987 

1983.986 

19BA.089 

1984.230 

1984.396 

1984.583 

1984.787 

1984.962 

1996.500 

13 

1979.000 

1979.000 

1979.000 

1979.000 

1983.857 

1983.958 

1984.108 

1984.273 

1984.455 

1984. 651 

1984.850 

1996.500 

lA 

1S78.900 

1973.900 

1978.900 

1978.900 

198 3.721 

1983.782 

1983.935 

1984.102 

1984.284 

1984.481 

1984. 690 

1984.701 

15 

1976.800 

1978.800 

1978.800 

1978.800 

1983.550 

1983. 575 

1983.719 

1983.885 

1984.070 

1984.272 

1984.490 

1984.530 

16 

1976. AOO 

1978. AOO 

1978. AOO 

1978. AOO 

1978. AOO 

1983.303 

1983.429 

1983.602 

1983.797 

1984.016 

1984.2 54 

1984.307 

n 

1978.000 

1978.000 

1978.000 

1978,000 

1978.000 

1982. 9A9 

1983.055 

1983.243 

1983.453 

1983.687 

1983. 94L 

1963.996 

IB 

1976.900 

1976.900 

1976.900 

1976.900 

1976.900 

1982.A96 

1982.584 

1982.797 

1983.032 

1983.284 

1983.555 

1903.614 

IS 

1975.800 

1975.800 

1975.800 

1975.800 

1975,800 

1981.982 

1982.043 

1982.292 

1982.550 

1982.817 

1983.099 

1983.149 

2C 

197A. 850 

197A.850 

197A.850 

197A.850 

197A.850 

1931.438 

1981.488 

1981.760 

1982.034 

1962.304 

1982.567 

1982.581 

21 

1973.900 

1973.900 

1973.900 

1973.900 

1973.900 

1980.929 

1980.950 

1981.241 

1981.535 

1981.834 

1982.136 

1993.500 

22 

1973.050 

1973.050 

1973.050 

1973.050 

1973.050 

1973.050 

1980.429 

1980.716 

1981.030 

1981.357 

1981.670 

1992.500 

23 

1972.200 

1972.200 

1972.200 

1972.200 

1972.200 

1972.200 

1979.888 

1980.155 

1980. a87 

1980.833 

1981.124 

1991.500 

2A 

1971,700 

1971.700 

1971.700 

1971.700 

1971.700 

1971.700 

1979.302 

1979.557 

1979.692 

1980.242 

1980.490 

1990.050 

25 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1978.667 

1978.914 

1979.242 

1979.594 

1979.708 

1 908.600 

26 

1971.200 

1971. 20U 

1971.200 

1971*200 

1971.200. 

1971.200 

1977.982 

1970.203 

1978.507 

1978.641 

1978.942 

1906.900 

27 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1971.200 

1977.162 

1977.379 

1977.672 

1977.997 

1978.038 

1905.200 

28 

1970.900 

1970.900 

1970.900 

1970.900 

1970.900 

1970,900 

1976.185 

1976.426 

1976.713 

1976.996 

1903.350 

1983.350 

2S 

1970.600 

1970.600 

1970.600 

1970. 600 

1970.600 

1975. 0A6 

1975.073 

1975.335 

1975.620 

1975.806 

1981.500 

1901.500 

30 

1969.650 

1959.650 

1969.650 

1969.650 

1969.650 

1973. 8A7 

1973.920 

1974.161 

1974.414 

1974.507 

1979,550 

1979.550 

31 

1968.700 

1968. 700 

1968.700 

1968.700 

1968,700 

1972.669 

1972.769 

1972.944 

1973.033 

1977.600 

1977.600 

1977.600 

32 

1968.350 

1968.350 

1968.350 

1968.350 

1960.350 

1971. 7A3 

1971.876 

1972.077 

1972.231 

1975.800 

1975.800 

1975.800 

33 

1968.000 

1968.000 

1968.000 

1968.000 

196 8.000 

1970.830 

1970.986 

1971.167 

1971.229 

1974.000 

1974.000 

1974.000 

3A 

1967.250 

1967.250 

1967,250 

196 7.250 

1969.801 

1969.892 

1970.064 

1970.224 

1972.500 

1972.500 

1972.500 

1972.500 

35 

1966.500 

1956.500 

1966.500 

196 6.500 

1968. 9A1 

1969.008 

1969.173 

1969.257 

1970.500 

1970.500 

1970.500 

1970.500 

36 

1965.750 

1965.750 

1965.750 

1965.750 

1968,093 

1968.212 

1960.431 

1968.442 

1969.000 

1969.000 

1969.000 

1969.000 

37 

1965.000 

1965.000 

1965.000 

1965.000 

1967.08A 

1967.235 

1967.312 

1967.500 

1967.500 

1967.500 

1967.500 

1967.500 

38 

196A.250 

196A. 2 50 

196A.2S0 

1966.113 

1966. 12A 

1966. 1A2 

1966.177 

1965.850 

1965,850 

1965.850 

1965.850 

1965.050 

39 

1963.500 

1963.500 

1953.500 

1965.833 

1965; 887 

1965.883 

1964.200 

1964.200 

1964.200 

1964.200 

1964.200 

1964.200 

AO 

1963.000 

1963.000 

1963.000 

1982.000 

1965. 767 

1982.000 

1963.000 

1963.000 

1963.000 

1963.000 

1963.000 

1963.000 
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ISGTHERMy^L LINE LOCATIONS 


I 

J 

T 

T-LOW 

T-HIGH 

FRAC 

2 

4 

1979 

1978.7146 

1979. 0198 

0.9352 

2 

6 

1980 

1979.6335 

1980.6155 

0.3732 

2 

7 

1981 

1980.6155 

1981.2554 

0,6009 

3 

6 

1981 

1 980.9097 

1931.4275 

0,1744 

3 

8 

1982 

1981-9712 

1982.0547 

0.3450 

4 

2 

1981 

1980,9280 

1981.0215 

0.7702 

4 

6 

19 82 

1981.7742 

1982.1509 

0.5995 

5 

4 

19 82 

1981, 876C 

1982.1379 

0,4734 

5 

7 

1983 

1982.7734 

1983.1401 

0,6178 

6 

6 

1903 

1982.9695 

1983.2837 

0,0971 

6 

8 

1984 

1983.6292 

1934. 0066 

0,9325 

7 

4 

1983 

1982.9004 

1933.1318 

0.4304 

7 

7 

1934 

1983.6831 

1934.0093 

0,9716 

8 

7 

19 84 

1933.9622 

1984.2495 

0.1317 

9 

6 

1984 

1983-9160 

1934.1392 

0.3764 

9 

10 

1985 

1984.9402 

1985-0405 

0,5961 

10 

5 

19 84 

1983.0633 

1984. 0369 

0.7876 

10 

10 

1985 

1984.9324 

1985,0713 

0.486B 

11 

5 

1934 

1933.9333 

1984.0806 

0.4292 

11 

10 

1985 

1984.8757 

1935. 0435 

0. 7409 

1.2 

5 

1984 

1983.9861 

1984.0894 

0.1348 

13 

6 

19 34 

1903.9583 

1934. 1084 

0.2 730 

14 

7 

1934 

1933.9353 

1984.1013 

0.3886 

15 

8 

1904 

1983.8845 

1934.0696 

0,6240 

16 

9 

1984 

1983.7971 

1984,0156 

0.9285 

17 

6 

1933 

1982.9492 

1983.0552 

0.4793 

18 

0 

1963 

1982.7974 

1933.0322 

0,S62B 

19 

6 

1982 

1981.9819 

1982.0435 

0.2537 

19 

10 

1983 

1902.8171 

1983,0991 

0.6485 

20 

8 

1982 

1981.7603 

1982. 0337 

0,3763 

21 

7 

19 81 

1900.9497 

1981.2412 

0,1725 

21 

10 

19 82 

1981.8337 

1982. 1360 

0,5501 

22 

8 

1981 

1980.7156 

1981, 0295 

0,9059 

23 

7 

1980 

1979-8884 

1980.1545 

0,4193 

23 

10 

1981 

1980. 8328 

1931,1240 

0.5742 

24 

9 

1980 

1979.8916 

1980, 2419 

0.3094 

25 

a 

1979 

1978.9136 

1979.2422 

0.2530 

26 

7 

1978 

1977.9819 

1978,2026 

0 o 0 8 1 9 

27 

10 

1978 

1977.9971 

1973.0378 

0-0719 

30 

7 

19 74 

1973.9197 

1974.1611 

0.3327 

31 

8 

1973 

1972.9436 

1973,0330 

0,6311 

32 

"7 

1972 

1971,8765 

1972.0774 

0.6148 

33 


1971 

1970.9861 

1971,1670 

0.0759 

34 

6 

1970 

1969.8923 

1970.0642 

0.6264 

35 

5 

1969 

1968.9409 

1969-0083 

0 « S 76 8 
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r; 


w > , 






MAX= 1998,50000 MIN= 1961.00000 
TEMP. INCREASES AS : 0-1^2-3-4-5-6-7-8-9-. 


pictorial view of temperature 


ui 

a\ 


0 




Xc 

Xc 

X: 


X: 


•»- 


y. 

9 

9 

9 

9 

9 

8 

0 


X-. 


Xs 


Xc 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

8 

0 




9 

9 

9 

9 

9 

9 

9 

9 

Q 

9 

9 

9 

Q 

✓ 

8 

8 

Q 


8 

8 

9 

9 

9 

9 

9 

9 

9 

9 

9 

a 

y 

9 

9 

a 

8 

3 

0 

8 

8 

8 

a 

9 

u 

y 

9 

9 

9 

9 

9 

9 

9 
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